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Dynamical  and  Structural  Studies  of  Molecules. 
NMR  Experiments,  Quantum  Chemstry,  and  Molecular 

Dynamics  Simulation. 


R.R.Emst,  T.  Bremi,  R  Brueschweiler,  S.Lienin,  R.Meyer,  P.Robyr, 

Ch.Scheurer,  andM.Utz 

Lab  orator  ium  fuer  Physikalische  Chemie, 

ETH-Zentrum,  8092  Zurich,  Switzerland 

NMR  experiments  provide  unique  insights  into  structural  features  and 
particularly  into  dynamical  aspects  of  molecules.  Chemical  shift,  J-coupling, 
and  relaxation  measurements  using  sophisticated  experimental  techniques 
form  the  basis  of  the  exploration.  The  analysis  of  the  data  requires  a  detailed 
knowledge  of  the  underlying  structural  dependence  of  the  NMR  parameters. 
Rarely  they  can  be  measured  with  sufficient  accuracy.  Often,  they  need  to  be 
calculated  by  quantum  chemical  methods.  The  interpretation  of  the 
intramolecular  dynamics  can  take  advantage  of  computer  modelling  by 
molecular  dynamics  simulations. 

Examples  that  combine  NMR  experiments,  quantum  chemical  calculations, 
and  molecular  dynamics  simulations  are  presented.  They  involve  intra¬ 
molecular  dynamics  of  peptides  and  proteins,  hydrogen  dynamics  in  transition 
metal  hydrides,  and  the  interpretation  of  local  order  in  amorphous  polymers. 
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Parallel  Calculation  of  NMR  Shieldings  on  Clusters  of 
Inexpensive  Computers:  Applications  to  Peptides 

Peter  Pulay 

Department  of  Chemistry  ad  Biochemistry 
Fulbright  College  of  Arts  ad  Sciences 
University  of  Arkansas,  Fayetteville,  Arkansas 

Strategies  for  efficient  parallelization  of  ab  initio  NMR  shielding  computations  at 
the  SCF  and  DFT  level  will  be  discussed,  with  emphasis  on  clusters  of 
inexpensive  computers.  Such  a  cluster  can  give  supercomputer  performance  at 
a  very  moderate  price,  opening  new  possibilities  for  the  computational 
approach.  Calculations  with  500-1000  basis  functions  can  be  routinely 
performed.  Large-scale  calculations  for  peptide  shielding  surfaces,  and  for 
polysilanes  will  be  described. 
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RECENT  DEVELOPMENTS  IN  THE  ACCURATE 
CALCULATION  OF  NMR  CHEMICAL  SHIFTS 


Jiirgen  Gauss 

Institut  fur  Physikalische  Chemie, 
Universitat  Mainz,  D-55099  Mainz,  Germany 


Recent  developments  for  the  accurate  calculation  of  NMR  chemical  shifts 
using  many-body  perturbation  theory  (MBPT)  and  coupled-cluster  (CC) 
methods  will  be  discussed.  In  particular,  we  will  present 

•  CC  calculations  for  a  series  of  small  diatomic  molecules  together 
with  a  treatment  of  rovibrational  corrections  and  temperature  ef¬ 
fects, 

•  a  chemical  application  (to  a  substituted  vinyl  cation)  where  only 
high-level  coupled-cluster  calculations  resolve  an  earlier  discrepan¬ 
cy  between  theory  and  experiment, 

•  a  recent  integral-direct  implementation  at  the  second-order  MBPT 
level  which  allows  electron-correlated  NMR  chemical  shift  calcula¬ 
tions  on  large  molecules  with  more  than  50  atoms. 
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THE  COUPLED-CLUSTER  TREATMENT  OF  NMR  AND  ESR  SPECTRA 
Rodney  J.  Bartlett,  Ajith  Perera,  Janet  Del  Bene,  Tom  Henderson 

Quantum  Theory  Project 
University  of  Florida 
Gainesville,  FL,  USA  32611 

Hydrogen  bonds  are  ubiquitous  in  chemistry  and  biochemistry,  being  essen¬ 
tial  to  life  itself.  However,  the  concept  of  a  hydrogen  bond  in  biochemistry 
consists  solely  of  a  double  well  potential,  regardless  of  the  atoms  and  circum¬ 
stances  involved  in  the  bond.  Consequently,  NMR  spectra  for  biomolecules  are 
invariably  interpreted  subject  to  this  underlying  model  of  an  H-bond.  We  will 
present  studies  of  NMR  chemical  shifts  and  coupling  constants  for  a  variety  of 
H-bonds  demonstrating  that  different  H-bonds  have  characteristic  NMR  spectra 
that  can  be  used  for  their  identification  and  detection.  We  also  consider  coupling 
constants  in  polypeptides,  presenting  purely  ab  initio  Karplus-type  relations  for 
the  conformations  of  peptide  units. 

In  ESR  we  consider  the  determination  of  hyperfine  coupling  constants  for 
quinones,  using  the  relaxed  density-matrix  approach  of  coupled-cluster  theory. 

*This  work  is  supported  by  the  US  Air  Force  Office  of  Scientific  Research 
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The  MCSCF  response  approach  to  relativistic  corrections  to 

NMR  parameters 

Kenneth  Ruud1 ,  Juha  Vaara2,  Olav  Vahtras3,  Boris  F.  Minaev3,  Hans  Agren3, 

and  Jukka  Jokisaari2 

1  Department  of  Chemistry,  University  of  Oslo,  P.O.Box-1033  Blindern, 

N-0315  Oslo,  Norway 

2  NMR  Research  Group,  Department  of  Physical  Sciences,  University  of 

Oulu,  P.O.Box  333,  FIN-9057 1  Oulu,  Finland 
3  Department  of  Physics  and  Measurement  Technology,  Linkoping 
University,  S-58183  Linkoping,  Sweden 


The  talk  will  outline  the  recent  progress  that  has  been  made  in  the  calcula¬ 
tion  of  relativistic  corrections  to  nuclear  magnetic  resonance  parameters  using 
linear  and  quadratic  MCSCF  response  theory.  The  talk  will  briefly  discuss  the 
theory,  with  emphasis  on  the  computational  strengths  and  weaknesses. 

The  main  part  of  the  talk  will  be  devoted  to  presenting  results  of  recent 
MCSCF  calculations  of  the  spin-orbit  corrections  to  the  nuclear  shielding  and 
indirect  spin-spin  coupling  constants  for  different  molecules.  Some  of  these 
studies  have  utilized  large  basis  sets  and  large  active  spaces,  providing  bench¬ 
mark  results  for  small  molecules  which  may  be  used  for  calibrating  more  ap¬ 
proximate  treatments  of  the  relativistic  effects. 

The  intemuclear  distance  dependence  of  the  nuclear  shieldings  and  the  spin- 
orbit  corrections  to  the  nuclear  shieldings  will  also  be  discussed. 
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ON  CHF  CALCULATIONS  OF 
SECOND-ORDER  MAGNETIC 
PROPERTIES  USING  THE  METHOD 
OF  CONTINUOUS 

TRANSFORMATION  OF  ORIGIN  OF 
THE  CURRENT  DENSITY 


S.  Coriani.,  P.  Lazzeretti,  M.  Malagoli  and  R.  Zanasi, 
Dipartimento  di  Chimica,  Universita  degli  Studi, 
via  G.  Campi  183,  41100  Modena,  Italy. 


Abstract 

A  fully  analytical  formulation  is  outlined  for  computing  molecular 
magnetic  susceptibilities  and  nuclear  magnetic  shieldings  via  a  contin¬ 
uous  change  of  origin  of  the  electronic  current  density  induced  by  an 
external  magnetic  field.  The  change  of  origin  is  described  in  terms  of 
a  (continuous)  arbitrary  shift  function  d(r).  Coupled  Hartree-Fock 
second-order  magnetic  properties  of  CH4  and  CO2  molecules  have 
been  computed,  using  the  special  choice  d(r)=  r  as  generating  func¬ 
tion.  A  detailed  analysis  of  results  obtained  with  a  variety  of  basis 
sets  reveals  that  such  a  method  is  not  as  good  as  previously  suggested. 
Large  basis  sets  must  be  used  to  obtain  accurate  magnetic  properties. 
On  the  other  hand,  all  the  components  of  theoretical  nuclear  magnetic 
shielding  evaluated  via  this  approach  are  independent  of  the  origin  of 
the  vector  potential.  In  general,  theoretical  magnetic  susceptibilities 
depend  linearly  on  the  distance  between  different  coordinate  frames, 
but  are  origin  independent  for  centre-symmetric  molecules. 
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A  Density  Functional  Theory  based  scheme  for  studying  ESR  and  NMR  properties  of 
small  molecules  embedded  in  larger  systems. 

Jacques  Weber  and  Tomasz  Adam  Wesolowski 
Departement  de  Chimie  Physique,  Universite  de  Geneve, 

30,  quai  Ernest- Ansermet,  CH-1211  Geneve  4,  Switzerland 

A  theoretical  formalism  for  studying  properties  of  molecules  embedded  in  larger 
systems  is  presented.  This  method  (Kohn-Sham  Equations  with  Constrained  Electron 
Density  [1-3],  KSCED)  is  based  on  a  modification  of  the  standard  method  of  Kohn  and 
Sham.  As  a  consequence  of  partitioning  the  total  electron  density  into  two  components 
(p=  pi+p2),  an  additional  functional  is  used  to  approximate  a  small  fraction  of  the  total 
kinetic  energy  (Tsnad  [pi,  p2]=  Ts[pi+p2]  -  Ts[pi]  -Ts[p2]). 

We  started  combining  KSCED  with  approaches  developed  recently  by  Malkin  et  al  [4] 
for  calculating  NMR  and  ESR  properties  using  density  functional  theory.  Our  present  paper 
will  presents  the  first  results  obtained  using  this  new  methodology: 
a)  the  influence  of  the  Ar  or  Ne  matrix  on  the  hyperfine  coupling  constant  of  Mg+,  b)  the 
effect  of  interactions  with  different  cations  on  the  C02"  hyperfine  coupling  constants.  The 
calculated  properties  are  in  excellent  agreement  with  experimental  results. 

[1]  T. A. Wesolowski  and  AWarshel  (1993)  J.  Phys.  Chem.,  97,  8050 

[2]  T. A. Wesolowski  and  J.Weber  (1996)  Chem.Phys.Lett,  248,  71 

[3]  T. A. Wesolowski  (1997)  J.  Chem.  Phys.,  106,  8516 

[4]  V.  Malkin,  O.L.  Malkina,  L.A.  Eriksson,  and  D.R.  Salahub,  in  “Theoretical  and 
Computational  Chemistry”  vol.  1,  “Density  Functional  Calculations”  P.  Politzer  and  J.M. 
Seminario,  Eds.,  Elsevier,  1995 
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NEW  QUANTUM  MECHANICAL  OPERATORS 
FOR  EPR  HYPERFINE  COUPLING  AND  NMR 
SPIN-SPIN  COUPLING  CONSTANTS 


Daniel  M.  Chipman 

Radiation  Laboratory,  University  of  Notre  Dame 
Notre  Dame,  IN  46556  USA 

and 

Vitaly  A.  Rassolov 

Department  of  Chemistry,  Northwestern  University 
Evanston,  IL  60208  USA 


New  global  operators  have  been  derived  to  express  the  Fermi  contact 
interaction  between  an  electron  and  a  nucleus  or  between  two  electrons.  These 
supply  a  useful  alternative  to  the  traditional  representation  in  terms  of  a  delta 
function  operator.  Among  their  interesting  formal  properties  are  proper  cusp 
behavior  of  the  density  even  when  the  underlying  wave  function  is  cuspless, 
such  as  is  the  case  with  Gaussian  basis  sets.  Certain  serious  difficulties 
known  for  the  behavior  of  the  global  operator  previously  introduced  by 
Hiller,  Sucher,  and  Feinberg,  such  as  improper  long  range  asymptotic 
behavior  and  very  inefficient  implementation  due  to  contributions  from  two- 
electron  integrals,  are  ameliorated  in  the  new  formulations. 

Practical  results  using  the  new  operators  are  presented  for  difficult  properties 
such  as  first-order  hyperfine  coupling  constants,  as  observed  in  EPR,  and 
second-order  indirect  nuclear  spin-spin  coupling  constants,  as  observed  in 
NMR.  Compared  with  the  traditional  delta  function  formulations,  insignificant 
additional  computer  time  is  required  to  carry  out  these  calculations.  It  is  found 
that  the  new  operators  can  give  higher  accuracy  from  medium  sized  basis  sets, 
and  also  faster  convergence  to  the  complete  basis  set  limit. 
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Characterization  of  Nuclear  Magnetic  Shielding, 
Indirect  Spin-Spin  Coupling  and  Electric  Field  Gradient 
Tensors  using  Solid-State  NMR. 


Rod  Wasylishen,  G.  Bernard ,  D.  Bryce,  K.  Eichele,  M.  Gee, 

M.D.  Lumsden,  S.  Kroeker,  G.  Ossenkamp,  R.  Schurko,  G.  Wu. 

Department  of  Chemistry,  Dalhousie  University,  Halifax,  Nova  Scotia, 

Canada,  B3H  4J3. 


Solid-state  NMR  is  ideally  suited  for  probing  the  orientation  dependence  of 
magnetic  shielding,  spin-spin  coupling  and  EFG  tensors.  NMR  studies  of 
single  crystals  provide  the  most  definitive  information  about  these  NMR 
tensors.  In  the  case  of  powder  samples,  it  is  advantageous  to  examine 
systems  where  the  crystal  symmetry  guarantees  that  the  NMR  parameters 
will  be  axially  symmetric.  Several  examples  demonstrating  the  utility  of 
solid-state  NMR  will  be  provided.  For  example,  nitrogen  shielding  tensors 
of  aniline-15N  and  pyridine-15N  cobaloximes;  phosphorus  shielding  tensors 
for  cis-  and  trans-substituted  diphosphetes,  Cu,  P  spin-spin  coupling  tensors 
for  Cu(I)-phosphine  complexes,  aluminium  shielding  and  EFG  tenors  for 
several  nearly  octahedral  complexes,  etc.  The  importance  of  first  principle 
calculations  will  be  discussed  and  several  challenging  problems  will  be 
presented. 
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SUBSTITUENT  EFFECTS  ON  TRANSITION-METAL 
CHEMICAL  SHIFTS. 

M.  Buhl 

Organisch-Chemisches  Institut,  Universitat  Zurich,  Winterthur erstr.  190, 

8057  Zurich,  Switzerland 

Chemical  shifts  of  transition-metal  nuclei  in  a  large  number  of  bonding 
situations  can  be  well  reproduced  with  modem  tools  of  density-functional 
theory.1  Effects  of  substituents  in  the  first  and  second  coordination  sphere 
can  be  rationalized,  e.g.  in  terms  of  geometrical  and  electronic  effects,  and 
the  parameters  governing  the  observed  trends  can  be  identified.  For  rhodium 
complexes  bearing  chelate  ligands,  for  example,  it  is  shown  that  the  Rh- 
ligand  distance  has  a  much  larger  effect  on  8(103Rh)  than  the  chelate  bite 
angle,  hitherto  assumed  to  be  of  prime  importance.2  The  computed 
coordination  geometry  and  51V  chemical  shift  of  the  vanadate-glycylglycine 
complex  below,  modeling  biologically  relevant  vanadate-peptide  inter¬ 
actions,3  is  indicated  to  be  quite  sensitive  to  the  pH  value. 


8(51V)calc  -502  ppm  H0  ^  -602  ppm 

1.  M.  Buhl,  M.  Kaupp,  V.  G.  Malkin,  O.  L.  Malkina,  J.  Comput.  Chern., 
submitted. 

2.  e.g.  J.  G.  Donkervoort,  M.  Buhl,  C.  J.  Elsevier,  Eur.  J.  Inorg.  Chem., 
submitted. 

3.  H.  Sigel,  A.  Sigel,  Eds.,  Metal  Ions  in  Biological  Systems,  Vol.  31, 
Marcel  Dekker,  New  York,  1995. 
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Calculations  of  Phosphorus  Nuclear  Magnetic  Shielding  Tensors 

Ulrich  Fleischer 

Department  of  Chemistry  and  Biochemistry,  University  of  Arkansas, 
Fayetteville,  Arkansas,  72701,  USA  -  Lehrstuhl  fiir  Theoretische  Chemie, 
Ruhr-Universitat  Bochum,  D-44780  Bochum,  FR  Germany 

The  nuclear  magnetic  shielding  of  a  nucleus  is  a  sensitive  probe  of  its  local 
environment.  In  addition  to  the  isotropic  value  the  nuclear  magnetic 
shielding  tensors  contain  some  spatially  resolved  information  and  are  thus 
even  more  sensitive  probes.  We  discuss  results  of  our  shielding  tensor 
calculations  on  phosphorus  compounds.  Phosphorus  has  not  only  favourable 
NMR  properties,  it  is  also  able  to  form  bonds  in  a  variety  of  different 
situations,  from  lowvalent  multiple  bonded  up  to  hypervalent  species. 
Electron  correlation  effects  are  more  pronounced  for  the  lowvalent  systems. 
HF-,  DFT-  and  MP2-data  of  a  few  systems  are  presented  and  discussed. 
Some  of  the  calculational  approaches  allow  for  a  breakdown  of  the 
calculated  tensors  into  contributions  of  localised  orbitals.  The  different 
orientations  of  the  principal  axes  system  and  the  variations  of  the  principal 
values  on  substitution  are  discussed  in  terms  of  interactions  of  occupied  and 
virtual  orbitals  induced  by  the  external  magnetic  field. 

To  study  the  influence  of  intermolecular  interactions  on  the  shielding 
tensors  calculations  on  clusters  were  performed.  The  structures  were  chosen 
to  model  the  local  environment  of  a  molecule  in  the  crystal.  The  results  are 
discussed  and  compared  to  solid  state  NMR  data. 
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Recent  Advances  in  the  Experimental  and 
Computational  Characterisation  of  Carbocations. 

Hans-Ullrich  SIEHL 

University  of  ULM,  Division  of  Organic  Chemistry  I 
D-89069  ULM,  Germany 

e-mail:  ullrich.siehl@chemie.uni-ulm.de,  FAX  +49-(0)73 1-502  2800 

Carbocation  chemistry  has  always  been  a  cornerstone  of  physical  organic  chemistry  and 
a  test  field  for  new  experimental  methods.  Recently  quantum  chemical  methods  for  the 
calculations  of  experimentally  observable  parameters  such  as  chemical  shifts  and  spin- 
spin  coupling  constants  have  become  available.  This  has  been  a  major  breakthrough  for 
the  structure  determination  of  reactive  intermediates.  We  have  applied  the  combination 
of  contemporary  experimental  high  field  NMR  methods  and  ab  initio  calculations  to  the 
investigation  of  new  carbocations  which  shed  some  new  light  on  controversial  problems 
and  demonstrate  the  fruitful  interplay  between  experiment  and  theory. 

Among  others  we  will  present  H,C  Cosy  NMR  spectra  of  silyl-substituted  hypercoordi- 
nated  bicyclobutonium  cations  and  accompanying  quantum  chemical  calculations  of 
geometries  and  NMR  chemical  shifts  and  coupling  constants. 


Si-Me 
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Theoretically  predicted  13C  NMR  Chemical  Shifts 
and  C-C  Coupling  Constants  as  Tools  for  the  Conformational  Analysis 
of  Natural  Products.  A  DFT  and  Molecular  Mechanics  Study 

Gemot  Frenking 

Fachbereich  Chemie,  Philipps-Universitat  Marburg, 
Hans-Meerwein-Strasse,  D-35032  Marburg,  Germany 


The  talk  will  focus  on  the  scope  and  limits  of  a  theoretical  approach  using 
Boltzmann  weighted  low-energy  conformations  given  by  molecular 
mechanics  in  conjunction  with  SOS-DFPT/IGLO  calculated  13C  NMR 

•j 

chemical  shifts  and  Jcc  coupling  constants  for  the  conformational  analysis 
of  polyketide  natural  products. 
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The  Challenge  of  EPR  Parameters 


Emil  Roduner 

Institute  of  Physical  Chemistry,  University  of  Stuttgart, 
Pfaffenwaldring  55,  D-70569  Stuttgart 

The  predictability  of  hyperfine  coupling  constants  of  small  organic  free 
radicals  has  reached  a  level  where  quantum  chemical  calculations  are  very 
useful  in  assisting  assignments  of  experimental  values  to  chemical 
structures.  However,  running  times  become  quickly  prohibitive  for  a 
meaningful  analysis  of  larger  systems,  or  when  temperature  dependencies 
require  extensive  averaging  over  vibrational  wave  functions. 

Transition  metals  play  an  important  role  in  the  chemistry  of  molecular 
complexes,  in  processes  essential  to  life,  and  in  catalysis.  Since  the  metal 
centres  are  often  paramagnetic,  EPR  becomes  the  method  of  choice  for  their 
investigation,  g- tensor  anisotropies  are  highly  sensitive  to  the  environment, 
hyperfine  anisotropies  and  Fermi-contact  interactions  reflect  the  distribution 
of  spin  densities.  The  interpretation  of  these  parameters,  however,  is  far  too 
often  still  persued  on  an  empirical  or  intuitive  level.  Rigid  quantum 
chemical  methods  which  include  electron  correlation  and  relativistic  effects 
seem  to  be  on  the  horizon  and  are  badly  needed  to  assist  the  experimentalist. 

Experimental  results  on  temperature  dependent  hyperfine  couplings,  on  g- 
and  hyperfine  tensors  of  transition  metal  ions  and  small  clusters,  and  on  ions 
interacting  with  an  organic  radical  will  be  reported.  Their  structures  and 
location  in  a  zeolite  lattice  are  expected  to  have  tremendous  impact  on  their 
catalytic  activity  and  selectivity.  Theoretical  support  is  needed  towards  a 
detailed  understanding  of  the  observations. 
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Applicability  of  DFT  and  Hybrid  HF-DFT  Methods  for  Evaluating 
EPR  Parameters  in  Bioradicals 


Leif  A.  Eriksson 
Dept,  of  Quantum  Chemistry 
Uppsala  University 
Box  518,  751  20  Uppsala,  Sweden. 
Email:  leifaxel@kvac.uu.se 


There  is  now  growing  evidence  that  radicals  play  a  crucial  role  in  several  im¬ 
portant  enzymatic  reaction  mechanisms,  both  in  terms  of  ’’storage  units”  for 
enzyme  reactivity  as  well  as  being  directly  involved  in  reaction  catalysis  and/or 
control.  During  the  last  few  years,  the  development  of  hybrid  and  gradient 
corrected  density  functional  theory  (DFT),  coupled  with  new  hardware  tech¬ 
nologies,  have  made  accurate  computational  studies  of  these  types  of  systems 
a  feasible  complementary  tool  to  experiments. 

In  the  present  seminar,  the  free  radical  governed  enzymatic  reaction  mech¬ 
anism  of  Ribonucleotide  Reductase  (RNR)  is  outlined,  based  on  theoretical 
calculations  at  the  B3LYP  level.  In  the  deoxyribonucleotide  generating  pro¬ 
tein  the  radical  reactivity  is  kept  ’’stored”  at  a  diiron-linked  tyrosine  residue 
some  35  away  from  the  active  site.  We  outline  the  role  of  the  tyrosine  radical, 
the  internal  radical  transfer,  the  radical  catalyzed  RNA  to  DNA  conversion, 
and  an  alternative  reaction  mechanism  observed  upon  subjecting  the  protein 
to  a  modified  azido-substituted  substrate,  leading  to  destruction  of  the  en¬ 
zymatic  function  and  reactivity. 

Several  examples  will  also  be  presented  of  other  bioradicals,  and  how  the  EPR 
parameters  evaluated  at  various  DFT  levels  can  serve  as  an  aid  for  experi¬ 
mental  interpretation  and  design  of  new  models.  We  will  in  this  section  focus 
on  DNA  and  sugar  radicals  generated  upon  exposure  of  DNA  to  ionizing  ra¬ 
diation,  and  quinone  radicals  involved  in  photosynthesis  in  green  plants  and 
blue-green  algae.  Results  obtained  by  combining  DFT  and  molecular  dynamics 
simulations  will  also  be  presented. 
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Ab  initio  calculations  of  ESR  parameters 

B.  Engels,  Inst,  fur  Phys.  und  Theor.  Chemie 
University  of  Bonn,  53115  Bonn,  Germany 


To  obtain  accurate  ESR  parameters  an  accurate  description  of  electron  correlation  is 
essential,  making  configuration  interaction  (Cl)  Or  other  high  level  approaches  mandatory 
in  the  case  of  ab  initio  calculations. .  On  the  other  hand  during  the  past  few  years  it  has 
been  shown  that  density  functional  theory  (DFT)  gives  also  accurate  ESR  parameters  in 
many  cases,  so  that  DFT  became  a  serious  alternative  to  conventional  approaches.  In  the 
first  part  of  this  contribution  the  accuracy  of  DFT  is  compared  with  methods  which 
combines  Cl  methods  with  perturbational  corrections.  New  developments  in  the 
individually  selecting  MR-CI  method  will  be  outlined. 

The  second  part  will  be  devoted  to  the  influence  of  the  nuclear  dynamic  on  ESR 
parameters. 
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HYPERFINE  COUPLING  DETERMINATION  FOR  FREE  RADICAL 
INTERMEDIATES  OF  PHOTOSYNTHETIC  ELECTRON  TRANSFER 


P.J  O'Malley 

Department  of  Chemistry,  UMIST,  Manchester,  M60  1QD,  UK. 

Free  radicals  of  quinones,  chlorophyll  and  tyrosine  are  formed  in  the  initial 
electron  transfer  reaction  of  photosynthesis.  Experimental  determination  of 
hyperfine  couplings  using  EPR/ENDOR  is  difficult  due  to  the  broadened 
spectra  observed  for  the  immobilised  species  involved.  We  have  recently  had 
success  in  the  quantitative  prediction  of  isotropic  and  anisotropic  hyperfine 
couplings  using  hybrid  density  functional  calculations(l,2).  These  are  proving 
to  be  an  invaluable  aid  in  the  assignment  of  experimental  spectra.  Examples 
illustrating  the  application  of  the  B3LYP  density  functional  method  to  the 
calculation  of  spin  densities  and  hyperfine  couplings  for  semiquinone,  tyrosyl 
and  bacteriopheophytin  radicals  will  be  presented. 

REFERENCES 

1.  P.J  O'Malley,  Density  functional  calculated  13-C,  17-0  and  1-H  isotropic  and 
anisotropic  hyperfine  couplings  for  the  plastosemiquinone  anion  radical  involved  in 
oxygenic  photosynthesis,  J.  Amer.  Chem.  Soc.,  1998  ,  120,  5093 

2.  P.J  O'Malley,  Effect  of  hydrogen  bonding  on  the  spin  density  distribution  and 
hyperfine  couplings  of  the  p-benzosemiquinone  anion  radical  in  alcohol  solvents:  A 
hybrid  density  functional  study  J.  Phys.  Chem.  A  1997,  101,  9813 
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Paramagnetic  Centers  in  Amorphous  Insulating  Films 

P.M.  Lenahan 

The  Pennsylvania  State  University,  University  Park,  PA  16802 


Microelectronic  technology  is  dominated  by  metal/oxide/silicon  (MOS) 
integrated  circuitry.  Advances  in  MOS  technology  involve  a  continuous 
decrease  in  the  dimensions  of  the  individual  devices  which  make  up  these 
integrated  circuits.  With  device  oxides  of  less  than  5  nm  in  production 
today,  the  technology  approaching  fundamental  physical  limits.  Point 
defects  in  amorphous  insulating  films  and  at  amorphous  film/silicon 
interfaces  play  important  roles  in  these  physical  limitations.  These  defects 
already  play  important,  not  always  deleterious,  roles  in  the  operation  of 
present  day  devices.  These  point  defects  have  been  studied  with  electron 
paramagnetic  resonance  (EPR).  EPR  studies  in  these  systems  are  fairly 
difficult.  If  the  films  studied  are  technologically  useful,  the  defect  densities 
are  usually  fairly  low.  Since  the  films  are  quite  thin,  one  typically  deals 
with  small  volumes  of  low  defect  density  material.  Although  these 
difficulties  have  somewhat  hampered  work  in  this  area,  about  twenty 
paramagnetic  centers  have  been  identified  in  these  systems.  Combinations 
of  "electrical"  and  magnetic  resonance  measurements  have  identified  the 
"electrical"  properties  of  most  of  these  centers.  My  presentation  will 
review  this  work. 
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THEORY  OF  ELECTRIC  FIELD  EFFECTS  ON  HYPERFINE  INTERACTION 

Shashi  P.  Kama 

Air  Force  Research  Laboratory 
Space  Electronics  and  Protection  Branch 
3550  Aberdeen  Ave,  SE 
Kirtland AFB,  NM  87 117-5776 

The  hyperfme  splitting  in  the  electron  spin  resonance  (ESR)  spectrum  of  a 
paramagnetic  species  experiences  a  linear  shift  in  the  presence  of  an  external  electric  field, 
provided  the  magnetic  nucleus  is  located  at  a  noncentrosymmetric  site.  The  electric-field- 
induced  linear  shift  is  also  observed  in  the  hyperfme  splitting  of  nuclear  quadrupole 
resonance  (NQR)  spectrum  of  a  nucleus  located  at  a  noncentrosymmetric  site  in  a 
molecule  or  in  crystal  lattice.  Thus,  the  linear  electric  field  effect  on  the  ESR  and  NQR 
hyperfme  splitting  provide  valuable  information  about  the  local  atomic  arrangement 
around  a  magnetic  nucleus.  The  magnitude  of  the  shift  depends  on  the  strength  of 
interaction  between  the  electron  spin  (charge)  distribution  and  the  external  magnetic  field, 
thus  providing  information  about  the  spin  (charge)  distribution  in  the  system. 

The  higher-order  electric-field  effects  are  not  generally  observed  because  of  their 
decreasing  magnitude.  However,  these  effects  could  become  important  and  observable 
when  external  electric  field  is  due  to  high  power  laser. 

This  paper  will  discuss  a  perturbation  theory  treatment  of  the  electric  field  effects 
on  ESR  and  NQR  hyperfme  couplings.  Theoretical  methods  to  calculate  the  electric  field 
effects  within  Hartree-Fock  approach  will  be  described.  In  order  to  illustrate  the 
applicability  and  relevance  of  the  theory,  a  few  examples  from  the  ab  initio  calculations  of 
the  first-order  effects  in  selected  molecules  and  radicals  will  be  presented. 
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A  GIAO-DFT  Approach  to  the  Calculation  of  NMR  and 
ESR  Parameters:  Quantitative  Results  and  Qualitative 

Interpretations: 

Georg  Schreckenbach,  Yosadara  Ruiz-Morales,  Stephen  K.  Wolff  and  Yang 
Kahndogin,  Sergei  Patchkovskii  and  Tom  Ziegler .  Department  of  Chemistry, 
University  of  Calgary, Calgary, Alberta,  Canada  T2N  1N4. 

We  present  new  implementations  and  applications  based  on  the  ADF 
(Amsterdam  Density  Functional  program  system  )  GIAO-DFT  scheme  for  the 
calculation  of  NMRla~c  and  ESR23  shifts  and  spin-spin  coupling  constants33. 
Results  will  be  presented  from  shift  calculations  based  on  a  new 
implementation  that  includes  spin  orbit-coupling ld  and  a  fully  relativistic 
treatment^  within  the  Zero-Order-Regular-Approximation4.  The  GIAO 
scheme  in  general  lends  itself  well  to  qualitative  interpretations.  This  will  be 
demonstrated  in  connection  with  calculations  on  organic,  inorganic  and 
organometallic  77Se4a  and  125Te4b  shifts  as  well  as  coordination  shifts  of 
hydride40,  CO44-0  and  phosphine4f  ligands  attached  to  a  metal  fragment. 

The  second  part  deals  with  the  calculation215  of  metal-ligand  NMR  spin 
coupling  constants  as  well  as  g-tensors  of  small  radicals.  It  will  be  shown 
that  environmental  effects  from  the  host  lattice  are  important  in  g-tensor 
calculations315  and  that  better  relativistic  methods  are  required  for  NMR  spin 
coupling  constants  involving  5d-elements.  An  account  of  our  NMR  and  ESR 
work  has  appeared6. 

References: 

1. (a)  G.Schreckenbach  and  T.Ziegler  J.Phys.Chem.  1995,99,606  (b)Schreckenbach, 
Georg,  and  Ziegler,  Tom.  Int.  J.  Quant.  Chem.  1996,69,753-766  (c)  G.Schreckenbach 
and  T.Ziegler  Int.  J.  Quant.  Chem.  1997,61,899.  (d)  Wolff, S.  and  Ziegler, T. 
J.Chem.Phys.  1998.  (e)  Wolff, S.  and  Ziegler, T.  JT.Chem.Phys.  sunmitted. 
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DENSITY  FUNCTIONAL  THEORY  CALCULATION 
OF  MAGNETIC  PROPERTIES:  IMPROVED 
FUNCTIONALS  AND  CURRENT  DEPENDENCE. 

S.M.  Colwell, 

Department  of  Theoretical  Chemistry,  University  of  Cambridge. 

In  order  to  calculate  magnetic  properties  of  molecules  using 
Density  Functional  Theory  it  is,  in  principle,  necessary  to 
incorporate  the  current.  Work  so  far  suggests  that  the  effect  of  the 
current  dependent  term  is  small,  and  that  it  is  more  important  to  use 
exchange-correlation  functionals  that  give  accurate  results  for  the 
orbital  eigenvalues. 

This  paper  will  survey  work  done  recently  in  Cambridge  to  test 
new  functionals  developed  by  Tozer  and  Handy  which  are  optimised 
for  these  eigenvalues.  It  will  also  consider  the  effect  of  the  current 
and  the  appropriate  form  of  the  current  dependent  functional  for 
molecular  densities. 
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Relativistic  calculations  of  the  NMR  shielding  and 


indirect  spin-spin  coupling  tensors 

Lucas  Visscher1,  Thomas  Enevoldsen2,  Trond  Saue3,  Hans-Jorgen  Aagard  Jensen2, 
Jens  Oddershede2  and  Gustavo  Aucar4 
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2.  Department  of  Chemistry,  Odense  University,  Campusvej  55,  DK-5230  Odense  M,  Denmark 

3.  IRSAMC  -LPQ,  Universite  Paul  Sabatier,  1 18,  route  de  Narbonne,  31062  Toulouse  CEDEX,  France 

4.  Physics  Department,  Exact  Science  Faculty;  UNNE,  Avda.  Libertad  5300,  Corrientes,  Argentina 


Abstract 

We  will  discuss  application  of  the  4-component  relativistic  RPA  method  to  include  relativistic 
effects  in  the  calculation  of  NMR  parameters. 

We  will  first  consider  the  interpretation  and  evaluation  of  the  DSO  contribution  in  a  relativistic 
formalism.  The  different  approximations  introduced  are  tested  in  calculations  of  the  spin-spin 
coupling  tensors  of  the  water  analogues  (H20,  H2S,  H2Se  and  H2Te). 

In  a  second  application  we  will  compare  benchmark  values  for  the  NMR  shieldings  and  indirect 
spin-spin  coupling  tensors  of  the  hydrogen  halides  with  the  result  of  other  approaches.  Here  we 
find  that  inclusion  of  the  spin-orbit  coupling  in  a  non-relativistic  formalism  suffices  to  recover 
the  relativistic  effect  on  the  hydogen  shieldings  but  fails  to  reproduce  the  much  larger  effect  on 
the  halogen  shieldings. 
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CALCULATIONS  OF  ELECTRONIC  G-SHIFTS 
BY  LINEAR  RESPONSE  THEORY  FOR  MULTI- 
CONFIGURATIONAL  SCF  WAVE  FUNCTIONS 


Maria  Engstrom,  Boris  Minaev,  Olav  Vahtras.  and  Hans  Agren 


Multiconfiguration  linear  response  theory  is  applied  to  the  calculation 
of  electronic  g-shifts,  normally  dominated  by  the  second-order  energy 
correction  involving  spin-orbit  and  orbital  Zeeman  effects.  Response 
theory  provides  an  alternative  method  for  evaluating  perturbational 
expressions,  where  summation  over  all  excited  states  is  achieved  im¬ 
plicitly  by  solving  linear  systems  of  equations.  Other  contributions 
(diamagnetic  and  relativistic)  are  first-order  corrections  and  are  eval¬ 
uated  as  expectation  values.  It  is  demonstrated  for  a  sequence  of  di¬ 
atomic  molecules  with  a  triplet  ground  state  what  experimental  agree¬ 
ment  can  be  expected  via  the  related  property  given  by  Curl’s  relation, 
the  spin-rotation  constants.  Improvements  of  the  present  formalism 
will  be  discussed  in  terms  of  non-variational  corrections  and  London 
orbitals. 
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Density  functional  calculations  of  ESR  parameters 
in  the  zero  order  regular  approximation  for  relativis¬ 
tic  effects 

Erik  van  Lenthe.  Paul  E.S.  Wormer  and  Ad  van  der  Avoird 

A  method  has  been  developed  for  the  calculations  of  ESR  param¬ 
eters  (g-tensor,  hyperfine  interaction  tensor)  of  Kramers  doublet  open 
shell  molecules,  which  uses  the  spinor  of  the  unpaired  electron  of  the 
paramagnetic  molecule,  obtained  from  a  density  functional  calculation. 

Spin-orbit  coupling  is  taken  into  account  variationally  using  the  re¬ 
cently  developed  zeroth  order  regular  approximated  (ZORA)  relativistic 
equation,  available  in  the  Amsterdam  Density  Functional  program  pack¬ 
age.  Relativistic  regular  approximated  expressions  for  the  Zeeman  term 
and  hyperfine  interaction  were  derived  and  implemented  in  the  ADF 
program,  with  the  use  of  numerical  integration  for  the  calculation  of  the 
corresponding  matrix  elements. 

Since  this  method  includes  the  spin-orbit  coupling  and  other  relativis¬ 
tic  effects  in  the  SCF  calculation,  and  uses  first  order  perturbation  theory 
for  the  Zeeman  terms  in  the  Hamiltonian,  it  is  not  much  more  expensive 
than  an  ordinary  SCF  calculation  on  a  molecule  without  external  field, 
and  can  easily  be  used  for  larger  molecules. 


L25 


PERSPECTIVES  ON  AB  INITIO  g-TENSOR 
CALCULATIONS 

Gerald  H.  Lushington 

Ohio  Supercomputer  Center,  1224  Kinnear  Road,  Columbus  OH,  USA 

43212-1163 


ABSTRACT 

Significant  progress  has  been  made  in  recent  years  toward  developing  vi¬ 
able  ab  initio  techniques  for  accurate  computation  of  electronic  g-tensors  for 
gas  phase  molecular  radicals.  Multi  reference  configuration  interaction  tech¬ 
niques  have  been  found  to  yield  highly  accurate  results  for  small  gas  phase 
molecules,  consistently  reproducing  general  trends  while  also  getting  good 
quantitative  agreement  (less  the  20%  deviation  from  experiment  for  larger 
g-shifts).  The  main  drawback  of  this  approach  is  its  high  computational 
cost,  a  factor  which  becomes  prohibitive  with  increasing  size  and  complexity 
of  the  system.  Much  more  computationally  tenable  methods  include  ROHF 
and  DFT  level  descriptions,  which  sacrifice  some  computational  accuracy  but 
generally  still  reproduce  basic  trends  in  g-shifts  and  bear  some  quantitative 
resemblance  to  experimental  values.  Many  radical  systems  of  interest,  how¬ 
ever,  do  not  lend  themselves  well  to  gas  phase  analysis  and  thus  are  often 
studied  in  solution,  matrix  or  solid  state  environments  -  conditions  under 
which  the  valence  character  of  the  radical  may  be  appreciably  perturbed, 
leading  to  significant  differences  in  g-shifts  as  compared  to  gas  phase.  This 
project  surveys  methodologies  and  results  for  conventional  (gas  phase)  cal¬ 
culations  and  provides  commentary  on  possible  condensed  phase  effects  and 
the  methods  that  may  be  applied  to  their  study. 
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Electron  Spin  Resonance  g  Tensor  Calculations 


Dylan  Jayatilaka 

Department  of  Chemistry,  The  University  of  Western  Australia, 
Nedlands  6009,  Australia 


Until  recently,  little  attention  has  been  given  to  calculating  the  g  tensor,  a  basic 
molecular  property  for  radicals  which  defines  the  magnetic  moment  in  an  applied 
external  magnetic  field. 

In  this  talk  I  will  describe  a  new  way  to  calculate  the  g  tensor  as  a  derivative  of  the 
energy  with  respect  to  an  applied  magnetic  field  [1].  Interestingly,  the  method  is 
formally  based  on  second  derivatives  of  the  square  of  the  Zeeman  splitting  energy. 
However,  the  formal  expression  can  be  simplified  to  a  first  derivative,  and  then  a 
new  Hellmann-Feynman  formula  for  the  g  tensor  can  be  derived.  To  get  results, 
a  General  Hartree  Fock  (GHF)  method  is  proposed  and  implemented.  The  GHF 
method  can  deal  with  the  required  spin-orbit  interactions  in  a  non-perturbative 
manner. 

g  tensor  elements  are  calculated  for  CO+,  CN,  MgF,  NO2,  NF2,  COj,  O3  ,  OHj, 
HCO,  and  H2CO+.  In  most  cases,  the  agreement  is  within  or  close  to  the  expected 
experimental  accuracy  of  500  ppm  for  neon  matrix  data,  but  there  are  still  some 
significant  discrepancies.  The  results  are  in  better  agreement  with  configuration  in¬ 
teraction  (Cl)  type  calculations  than  density  functional  theory  (DFT)  calculations, 
although  the  differences  are  small.  The  results  are  significantly  different  from  per¬ 
turbation  theory  predictions.  It  is  found  that  the  diamagnetic  contributions  to  the 
g  tensor  are  not  significant  to  within  the  experimental  errors,  if  matrix  isolation 
data  is  used. 

The  g  tensor  for  C0CI4-  ion  is  also  studied,  to  demonstrate  that  the  method  can  be 
applied  to  non-doublet  systems. 

References 

1.  “Electron  Spin  Resonance  g  Tensors  from  General  Hartree  Fock  Calculations” 
D.  Jayatilaka,  J.  Chem.  Phys  108,  7587  (1998). 
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APPLICATION  OF  NUCLEAR  SHIELDING  SURFACES  TO  THE 
FUNDAMENTAL  UNDERSTANDING  OF  ADSORPTION  AND  DIFFUSION 

IN  MICROPOROUS  SOLIDS 

Cynthia  Jameson 
University  of  Illinois  at  Chicago 


The  intermolecular  chemical  shift  provides  a  means  by  which  very  detailed 
questions  can  be  answered  towards  a  fundamental  understanding  of  two  basic  processes  in 
the  many  technological  applications  of  microporous  materials:  adsorption  and  diffusion. 
We  will  demonstrate  that  the  detailed  equilibrium  properties  of  physisorbed  small 
molecules  and  the  rates  of  cage-to-cage  migration  in  zeolites  can  be  observed  and 
quantitatively  reproduced  by  theoretical  calculations.  Some  of  the  questions  which  have 
been  answered  include:  When  molecules  are  adsorbed  in  a  microporous  solid  at  a  given 
loading,  what  fraction  of  the  cavities  have  exactly  n  molecules?  Is  adsorption  of  one  type 
of  molecule  enhanced  or  diminished  by  competition  with  another  type?  How  many 
molecules  of  type  2  can  be  found  in  those  cavities  that  have  exactly  n  molecules  of  type  1? 
How  do  the  type,  size,  and  locations  of  cations  affect  the  distribution  of  molecules  within  a 
microporous  solid?  How  is  the  NMR  chemical  shift  affected  by  confinement  with  n 
molecules  of  another  type  in  a  cavity?  A  zeolite  cage  provides  a  model  system  for  the 
understanding  of  solvent  effects  on  chemical  shifts,  for  the  theoretical  investigation  of  the 
temperature  dependence  of  medium  chemical  shifts  under  constant  volume  conditions  in 
which  very  detailed  comparisons  with  experiment  may  be  carried  out. 
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NMR  views  of  Proton  Transfer  and  Hydrogen  Bonding 

Hans-Heinrich  Limbach 

Freie  Universitdt  Berlin,  Institut  fur  Organische  Chemie, 

Takustr.  3,  D-14195  Berlin,  Germany 

The  problem  of  localization  and  transfer  of  hydrogen  isotopes  in  hydrogen 
bonded  hetercyclic  compounds,  acid-base  complexes,  transition  metal  hydrides 
embedded  in  liquids  and  solids  is  a  challenging  problem  of  chemistry.  In  this 
contribution  liquid  and  solid  state  NMR  -  methods  are  presented  which  are 
especially  adapted  for  the  study  of  these  problems.  For  certain  model  systems 
kinetic  hydrogen/  deuteron/tritium-isotope  effects  on  proton  transfer  through  and 
across  a  barrier  can  be  obtained  by  various  methods  of  dynamic  NMR  spectroscopy, 
ranging  from  the  second  to  the  nanosecond  timescale.  In  the  case  of  low-barriers 
hydrogen  bonds  H/D-isotope  effects  on  the  hydron  locations  are  observed  using 
dipolar  solid  state  NMR  which  are  related  to  H/D  isotope  effects  on  chemical  shifts. 
Novel  low-temperature  NMR  methods  allow  us  to  study  acid-base  interactions  in 
polar  solvents.  Zero-point  energies  are  directly  measured  via  H/D  fractionation 
factor  measurements.  For  the  first  time,  scalar  couplings  between  heavy  atoms  of 
hydrogen  bonds  are  observed  which  indicate  the  covalent  character  of  these  bonds. 
All  experimental  values  are  used  to  model  theoretically  the  details  of  the  proton 
motion  in  hydrogen  bonds. 
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The  SOS-DFPT  method  implemented  in  the  deMon-NMR  program  has 
proven  to  be  an  efficient  and  highly  accurate  source  of  NMR  chemical 
shifts.  Progress  will  be  reported  on  applications  involving  complex 
systems. 

•  Hydrogen-bonded  imidazole  complexes  were  chosen  to  model  the 
catalytic  triad  of  serine  proteases.  The  focus  is  on  nitrogen  atoms  as  they 
are  known  to  be  extremely  sensitive  to  intermolecular  interactions  and  to  the 
protonation  state.  The  validity  of  various  aspects  of  the  models  can  be 
established  using  SOS-DFPT. 

•  29Si  and  27 A1  shifts  have  been  calculated  for  models  of  the  nine 
different  sites  in  zeolite- {3.  The  Si  shifts  for  single-T-site  models  already 
allow  an  assignment  with  a  relative  error  less  than  1  ppm  and  an  absolute 
error  of  about  5ppm.  Two-site  models  reduce  the  absolute  error  to  1  ppm. 
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ELECTRIC  FIELD  GRADIENTS  IN  SOLIDS  BY  THE 
FULL-POTENTIAL  LAPW  METHOD 

K.  Schwarz  and  P.BIaha 

Technische  Universitat  Wien,  A- 1060  Vienna,  Austria 

NMR  is  one  of  the  experimental  techniques  to  measure  electric  field  gradients 
(EFGs)  by  nuclear  quadrupole  interaction.  The  EFG  is  a  traceless  tensor 
defined  as  the  second  derivative  of  the  electrostatic  potential  with  respect  to 
coordinates.  Up  to  about  10  years  ago  EFG  experiments  were  often  interpreted 
with  a  simple  point  charge  model  from  which  the  electrostatic  potential  and 
thus  the  EFG  is  deduced,  but  Stemheimer  antishielding  factors  were  needed  to 
correct  for  the  shortcoming.  Modem  electronic  structure  calculations  provide 
among  other  quantities  the  electron  density  distribution  from  which  the  EFG 
can  be  calculated  without  further  approximations.  For  solids  such  calculations 
are  based  on  density  functional  theory  (DFT)  according  to  Kohn  and  Sham. 
The  interacting  system  is  mapped  in  a  unique  manner  onto  an  effective  non¬ 
interacting  system  with  the  same  total  density.  The  effective  local  one-particle 
potential  consists  of  a  mean-field  (Hartee)  part  and  an  exchange-correlation 
part  which,  in  principle,  incorporates  all  correlation  effects  exactly.  However, 
the  exact  functional  form  of  this  potential  is  not  known  and  thus  one  needs  to 
make  approximations.  The  local  density  approximation  (LDA)  is  often  used 
but  recently  an  improvement,  namely  the  generalized  gradient  approximation 
(GGA)  has  reached  (almost)  “chemical  accuracy”. 

One  among  the  most  accurate  schemes  for  solving  the  Kohn  Sham 
equations  is  the  full-potential  linearized  augmented  plane  wave  (LAPW) 
method.  In  this  method  the  unit  cell  is  partitioned  into  (non-overlapping) 
atomic  spheres  centered  around  the  atomic  sites  and  an  interstitial  region, 
where  in  the  former  atomic-like  functions  and  in  the  latter  plane  waves  are 
used.  In  FP-LAPW  the  potential  and  charge  density  are  expanded  into  lattice 
harmonics  (inside  each  atomic  sphere)  and  in  a  Fourier  series  (in  the 
interstitial  region)  and  thus  they  are  completely  general.  The  numerical  (radial) 
wave  functions  guarantee  a  very  accurate  description  of  the  electron  density 
distribution,  especially  near  the  nucleus  and  its  anisotropy  that  leads  to  the 
EFG.  Application  to  insulators,  metals,  magnets,  superconductors,  etc  lead  to 
insight  into  the  origin  of  EFGs  and  good  agreement  with  experiments. 
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hyperfine  interaction  in  solids 
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Institute  for  Physical  Chemistry,  University  of  Munich,  Theresienstr.  37, 
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A  linear  response  formalism  will  be  presented  that  has  been  developed  within 
the  framework  of  spin-polarized  relativistic  multiple-scattering  theory  (SPR- 
KKR)  and  the  Green’s  function  technique  to  allow  for  a  calculation  of  mag¬ 
netic  susceptibilities  and  Knight-shifts.  This  very  general  approach  allows 
to  deal  with  paramagnetic  as  well  as  ferromagnetic  solids.  In  addition  the 
approach  is  suited  to  deal  with  all  spin  and  orbital  contributions  to  the 
magnetic  susceptibility  and  Knight-shift  on  the  same  level.  Furthermore  all 
relativistic  effects  are  taken  properly  into  account.  As  a  consequence  one 
finds  -  apart  from  the  conventional  orbital  Van  Vleck-contribution  to  the 
susceptibility  and  Knight-shift  -  an  additional  orbital  contribution  due  to 
the  spin-orbit  interaction.  Applications  will  be  presented  for  the  pure  3d-, 
4d-  and  5d-transition  metals.  The  extension  of  the  presented  scheme  to  al¬ 
loys  and  inhomogeneous  systems  as  for  example  surfaces  is  straightforward 
and  will  be  discussed  in  short.  In  addition  an  analogous  fully  relativistic 
treatment  of  other  hyperfine  interaction  parameters  of  solids  as  for  example 
the  nuclear  spin-lattice  relaxation  time  will  be  sketched. 

The  second  part  of  the  contribution  will  be  devoted  to  spontaneously  mag¬ 
netized  solids.  For  the  corresponding  hyperfine  field  a  proper  decomposition 
of  these  fields  into  the  conventional  Fermi-contact,  spin-dipolar  and  orbital 
parts  will  be  presented  and  discussed.  In  discussing  corresponding  results, 
the  main  emphasis  will  be  laid  on  the  spin-orbit  induced  orbital  contribu¬ 
tions  to  the  hyperfine  field  that  are  by  no  means  neglectable  compared  to  the 
dominating  Fermi-contact  part.  Finally  results  of  calculations  of  the  spin- 
orbit  induced  electric  field  gradients  in  cubic  systems  will  be  discussed  to 
demonstrate  the  symmetry  breaking  effect  of  the  spin-orbit  interaction. 
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I  will  present  a  theory  for  the  ab-initio  computation  of  the  NMR  chem¬ 
ical  shifts  in  condensed  matter  sy stems  [1].  The  formalism  can  be  ap¬ 
plied  to  periodic  crystals  or,  using  a  super-cell  technique,  to  liquids  and 
amorphous  materials.  I  will  illustrate  the  method  with  applications  to 
crystalline  and  amorphous  silica,  to  amorphous  carbon  thin  films,  and  to 
ice  and  liquid  water. 

[1]  F.  Mauri,  B.  Pfrommer,  and  S.  G.  Louie,  Phys.  Rev.  Lett.  77, 
5300  (1996). 
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ISOLATED  WATER  MOLECULE. 

NEW  EXPERIMENTAL  AND  THEORETICAL  DATA 

N.  M.  Sergeyev 

NMR  Laboratory,  Department  of  Chemistry 
Moscow  State  University,  119899,  Moscow  Russia 

During  last  years  new  experimental  and  theoretical  results  were  obtained  for  NMR 
parameters  of  water  molecule.  Experimental  data  were  obtained  for  dilute  solutions 
of  water  in  nitromethane  where  proton  exchange  is  stopped  and  water  can  be 
postulated  to  exist  in  monomeric  form  or  as  isolated  molecules.  The  H-D  coupling 
constant  and  170-H  and  170-D  coupling  constants  in  water  were  obtained  with  high 
accuracy  of  about  0.01  Hz.  This  allowed  us  to  find  D/H  isotope  effects  on  coupling 
constants  and  D/H  induced  170  isotope  shifts.  The  primary  isotope  effect  on  170-H 
coupling  constant  was  found  to  negative  in  accord  with  theoretical  estimates.  Also  a 
small  but  significant  nonadditivity  of  D/H  isotope  effects  for  170  chemical  shifts  was 
established.  Dilute  solutions  of  water  in  nitromethane  also  lead  to  long  proton  and 
deuteron  relaxation  times.  For  mixtures  of  H20  and  D20  the  equilibrium  constant  for 
the  reaction  H20  +  D20  =  2  HDO  was  found  using  170  NMR  and  also  a  combination 
of  'H  and  2D  NMR  spectra.  Accuracy  for  experimental  and  theoretical  results  is 
analyzed  and  the  sense  of  the  comparison  between  experimental  and  calculated  is 
discussed.  At  present  in  some  cases  calculated  NMR  parameters  for  small  molecules 
can  be  considered  as  more  accurate  than  the  experimental  ones.  It  is  worth  noting 
that  the  isotope  effects  studies  are  an  important  step  in  checking  any  theoretical 
approach  as  the  isotope  effects  on  NMR  parameters  can  be  considered  as  small 
deviations  from  the  equilibrated  values  and  can  be  treated  in  terms  of  the  theory  of 
perturbation. 
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Isotope  effects  on  the  NMR  parameters  of  water  calculated  from 
correlated  nuclear  shielding  and  nuclear  spin-spin  coupling  surfaces 

W.T.  Raynes 

Department  of  Chemistry,  University  of  Sheffield,  Sheffield,  S3  7HF,  UK 


Coefficients  of  the  two  nuclear  shielding  surfaces  cr(O)  and  a(H)  and  the  two 
nuclear  spin-spin  coupling  surfaces  J(H,H)  and  J(0,H)  of  the  water  molecule 
have  been  calculated  to  second  order  in  internal  displacement  coordinates  using 
a  correlated  wavefunction.  These,  coefficients  have  been  used  to  calculate  the 
effects  of  isotopic  substitutions  (!H,  2H  and  160,  170,  180)  on  the  four  NMR 
parameters  over  a  range  of  temperature.  The  results  are  in  excellent  agreement 
with  recent  experimental  measurements  on  water  isotopomers  dissolved  in 
"inert"  solvents.  Analysis  of  the  calculated  results  shows  that  bond  stretching  at 
first  order  dominates  the  nuclear  motion  corrections,  but  that  both  bond 
stretching  and  angle  bending  at  second  order  make  substantial  contributions. 
The  observed  nonadditivity  of  deuterium  substitution  on  the  170  isotope  shift 
can  be  largely  attributed  to  second  order  bending.  Some  unexpected  features  of 
the  shielding  and  spin-spin  coupling  surfaces  of  several  molecules  will  also  be 
discussed. 
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QED  THEORY  FOR  MAGNETIC  RESPONSE  PROPERTIES  IN 

MOLECULES 

RODOLFO  H.  ROMERO  AND  GUSTAVO  A.  AUCAR 


Abstract.  The  current  way  to  include  relativistic  effects  on  calculations  of 
magnetic  molecular  properties  are  separated  in  scalar,  2-component  or  4- 
component  (full  relativistic)  based  formalism.  All  of  them  starts  from  Dirac, 
Pauli  or  Schrodinger  schemes  and  use  perturbation  theory  (or  response  the¬ 
ory)  to  include  the  correction  to  the  energy  which  is  related  with  the  property 
analyzed.  In  this  communication  we  present  a  new  formalism  which  treat 
quantized  magnetic  field  and  second  quantized  boson  and  fermion  operators; 
this  is  done,  to  our  knowledge  for  the  first  time.  We  use  perturbation  theory 
applied  also  on  operators.  The  fine  structure  constant  is  used  as  perturbative 
parameter.  This  formalism  is  applied  on  NMR-J  tensor  calculations.  We  ob¬ 
tain  the  different  contributions  to  J,  i.e.,  the  usual  direct  and  indirect  nuclear 
spin  couplings  when  magnetic  fields  are  assumed  as  classical.  We  also  get  the 
QED  contributions  in  a  natural  manner.  Relativistic  corrections  arise  when 
their  4-component  implementation  is  reduced  by  current  methods  to  two-  or 
one-component  expressions. 


Facultad  de  Ciencias  Exactas  y  Naturales  y  Agrimensura.  Universidad  Nacional  del 
Nordeste.  Avda.  Libertad  5300  (3400)  Corrientes,  Argentina 
E-mail  address:  rhromeroGexa . unne . edu . ar ,  gaafirec.unne.edu.ar 
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CALCULATIONS  OF  NMR  PARAMETERS:  SOME 
APPLICATIONS  TO  LARGE  ORGANIC  SYSTEMS 

Jin  Czemek  and  Vladimir  Sklenar 
Laboratory  of  Biomolecular  Structure  and  Dynamics,  Faculty  of  Science, 
Masaryk  LIniversity,  Kotlarska  2,  CZ-61 1  37  Brno,  CZECH  REPUBLIC, 
e-mail  czernek@chemi.muni.cz,  sklenar@chemi.muni.cz 

For  large  (containing  16  nonhydrogen  atoms)  organic  molecules  -  2,3’-, 
2,5’-,  and  3’,5’-anhydrothymidines  -  chemical  shifts  were  obtained  using 
SOS-DFPT-IGLO1,  UDFT-GIAO2,  and  CPHF-GIAO3  methods.  Shielding 
for  both  experimental  (based  on  measured  coupling  constants  and 
subsequent  MM  optimizations)  and  RHF/6-31G**  geometries  was 
investigated.  Rather  good  agreement  with  experimental  values  of  13C  (rms 
deviations  of  relationships  between  theoretical  and  measured4  data  for  these 
methods  were  5.5,  5.2,  and  6.6  ppm,  respectively)  and  !H  chemical  shifts 
(0.233,  0.269,  and  0.297)  will  be  demonstrated.  SOS-DFPT-IGLO  method, 
being  the  most  reliable,  was  more  than  3  times  faster  than  G I  AO-based 
approaches. 

SOS-DFPT5  was  employed  to  calculate  selected  *H  -  *H  spin-spin 
coupling  constants  in  these  systems.  As  will  be  shown,  excellent  agreement 
with  experiment  can  be  obtained  which  facilitates  structural  studies  of 
anhydrothymidines4. 
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NMR  calculations  on  Hydrogen  Bonded  Imidazole 
Complexes-  Application  to  Serine  Proteases 

Michaela  Flock1,2  and  Hung  Guo3  and  Chung  Wong4  and  Dennis  Salahub2,3 
1  Laboratory  of  Quantum  Chemistry,  University  of  Leuven,  Celestijnenlaan  200 
P,  B-3001  Heverlee-Leuven,  BELGIUM 
2  De  partement  de  chimie,  Universite  de  Montre  al,  C.P.  6128,  succursale 
Centre-ville,  Montre  al,  Que  bee,  H3C  3J7,  CANADA 
3  CERCA-Centre  for  Research  on  Computation  and  its  Applications,  suite  400, 
5160  boul.  De  carie,  Montre  al,  Que  bee,  H3X  2H9,  CANADA 
4  Mount  Sinai  School  of  Medicine,  Department  of  Physiology  and  Biophysics, 

New  York,  NY,  10029-6574  USA 


The  imidazole  ring,  which  forms  the  side  chain  of  the  amino  acid  histidine, 
plays  an  important  role  for  protein  structure  stabilisation  and  is  involved  in  a 
number  of  enzymatic  reactions,  due  to  its  ability  to  form  hydrogen  bonds  via 
the  two  different  nitrogen  atoms  and  even  to  act  as  acid  or  base.  In  serine 
proteases  the  imidazole  ring  is  part  of  the  catalytic  triad,  which  consists  of  the 
via  hydrogen  bonds  connected  side  chains  of  aspartate,  histidine,  and  serine. 
Proton  and  nitrogen  NMR  techniques  where  applied  to  gain  information  about 
hydrogen  bonding  and  protonation  state  of  the  imidazole  ring,  in  order  to  support 
one  of  the  proposed  catalytic  mechanisms. 

We  studied  the  nitrogen  chemical  shifts  in  neutral  and  protonated  imidazole 
rings,  and  in  a  number  of  their  hydrogen  bonded  complexes.  Furthermore  we 
present  the  results  of  our  chemical  shift  calculations  on  the  active  site  of  serine 
proteases  using  models  of  increasing  size. 


P4 


CALCULATION  OF  DEUTERIUM  ISOTOPE  EFFECTS 
ON  13C  CHEMICAL  SHIFTS  IN  BENZOCYCLOBUTENE 

Z.  Glasovac.  P.  Novak,  Z.  Meic  and  M.  Eckert-Maksic 

Department  of  Chemistry,  Ruder  Boskovic  Institute,  HR- 10001  Zagreb,  CROATIA 


As  a  part  of  our  investigation  of  the  molecular  and  electronic  structure 
of  benzocyclobutene  (l)1  and  its  derivatives,  we  report  here  on  the  deuterium 
isotopic  effect  (DIE)  on  13C  NMR  chemical  shifts  in  2-5.  The  measured  DIE 
values  are  correlated  with  changes  of  net  charges  on  carbon  atoms  and  with 
changes  in  isotropic  magnetic  shielding  calculated  by  using  GIAO-B3LYP 
method  employing  several  basis  sets. 


1  X  =  X'  =  Y  =  Y'  =  H 

2  X  =  D;X  =  Y  =  Y  =  H 

3  X  =  X=D;Y  =  Y  =  H 

4  X  =  Y  =  D;X  =  Y'  =  H 

5  X=X=Y=Y=D 


Presence  of  the  deuterium  was  modeled  by  shortening  of  the  C-D  bond 
distance  by  0.012  A  relative  to  the  optimized  length  in  the  corresponding  C-H 
bond(s).  Theoretical  results  compare  favorably  with  the  experimental  data. 


1.  M.  Eckert-Maksic,  Z.  B.  Maksic,  and  M.  Klessinger,  J.  Chem.  Soc.,  Perkin  Trans.  2, 
(1994)  285;  Z.  B.  Maksic,  M.  Eckert-Maksic,  M.  Hodoscek,  W.  Koch,  and  D.  Kovacek, 
in:  Molecules  in  Natural  Science  and  Medicine,  An  Encomium  for  Linus  Pauling,  Z.  B. 
Maksic  and  M.  Eckert-Maksic  (Eds.)  Ellis  Horwood,  Chichester,  1991,  p.  333. 
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Prediction  of  29Si  MAS  NMR  chemical  shifts  in  zeolites 

P.  PULLUMBI(a) ,  A.  GOURSOT(b) ,  G.  VALERIO(b),  R.  VETRIVEL(b) 

and  D.  R.  SALAHUB(d) 

(a)  Air  Liquide,  Les  Loges  en  Josas,  Jouy  en  Josas,  France 

(b)  UMR  5618  CNRS,  Ecole  Nationale  Superieure  de  Chimie,  8  rue  de  l'Ecole  Normale 

34296  Montpellier  Cedex  5,  France 

(c)  Catalysis  Division,  National  Chemical  Laboratory,  411008  Pune,  India 

(d)  Departement  de  Chimie,  Universite  de  Montreal,  CP  6128,  Succursale  Centre  Ville, 

Montreal,  Quebec  H3C3J7,  Canada 


29Si  NMR  chemical  shifts  have  been  evaluated  from  the  NMR  shielding  tensors 
calculated  both  within  the  DFT  and  HF  formalisms.  The  Si  NMR  chemical  shifts 
have  been  calculated  for  the  nine  crystallographically  distinct  Si  sites  of  the  zeolite-p 
lattice.  The  calculations  were  carried  out  on  cluster  models  Si(OSiH3)4  and  R3Si- 
OSiR3  (R=OSiH3)  representing  one  (IT)  and  two  (2T)  sites  in  the  zeolite  lattice, 
respectively.  Using  the  IT  models,  the  nine  signals  of  the  experimental  spectrum  have 
been  assigned  with  a  relative  error  of  less  than  lppm,  the  average  absolute  error  being 
estimated  at  about  3-5  ppm.  The  use  of  a  "fragment-averaging"  technique  based  on  the 
results  obtained  with  the  2  T  models  leads  to  calculated  absolute  NMR  shifts  with  an 
accuracy  of  ±  1  ppm  for  DFT  and  ±  2  ppm  for  HF. 
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EXPERIMENTAL  AND  THEORETICAL  TENSORS  OF 
31P  AND  77Se  NUCLEAR  MAGNETIC  SHIELDING  IN 
(«£0Pe-O)2(S)PSeSeP(S)(O-«£0Pe)2 

G.  Grofmann.  D.  Scheller,  G.  Zahn,  Dresden  University  of  Technology, 

Germany 

O.L.  Malkina.  V.  G.  Malkin.  Slovak  Academy  of  Sciences,  Bratislava, 

Slovakia 

H.  Schmitt,  U.  Haeberlen,  Max-Planck-Institute  of  Medical  Research, 

Heidelberg,  Germany 

The  principal  values  of  P  and  Se  shielding  tensors  were  obtained  by 
Potrzebowski1  in  solid-state  MAS  NMR  investigations  of  the  compound 
mentioned  above.  The  orientation  of  the  principal  axes  of  the  tensors  with 
respect  to  a  molecular  coordinate  system  remains  hidden  by  investigations 
of  crystalline  powders.  This  missing  information  can  be  obtained  by  NMR 
goniometer  measurements  of  a  monocrystal  or  by  quantum  mechanical 
nuclear  shielding  calculations.  In  this  paper  we  will 

-  compare  the  complete  experimental  and  calculated  P  tensors, 

-  compare  the  experimental  and  calculated  principal  values  of  Se  tensors, 

-  try  to  understand  the  intensity  relations  of  two  NMR  signals  in  the  MAS 
and  goniometer  spectra. 

For  this  purpose,  3fP  goniometer  NMR  spectra  of  a  “monocrystal”  were 
measured  and  analyzed,  SOS-DFPT  IGLO  calculations  of  the  molecule 
mentioned  were  carried  out  and  the  crystal  structure  suggested  in  the 
literature2  as  a  triclinic  one  was  reinvestigated. 

The  same  conclusions  follow  from  the  monocrystal  NMR  investigation  as 
well  as  from  the  SOS-DFPT  IGLO  calculation:  the  principal  axis  3  deviates 
only  slightly  from  the  P=S  direction,  the  principal  axis  2  is  nearly 
perpendicular  to  the  P-Se  bond.  The  intensity  relations  of  the  NMR  signals 
can  be  understood  by  the  presence  of  a  monoclinic  phase  in  addition  to  the 
dominant  triclinic  phase. 


1  M.  J.  Potrzebowski.  Magn.  Resort.  Chem.  35,  8  (1995). 

2  M.  J.  Potrzebowski,  J.  Blaszczyk,  M.  W.  Wieczorek,  J.  Baraniak,  W.  J.  Stec,  in  press: 
Solid  State  NMR. 
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NMR  calculations  of  truly  large  systems  with  a  simplified 
DFT-LDA-Tight-Binding  scheme 

Th.  Heine1,2.  G.  Seifert1,  and  P.W.  Fowler2 
1  Institut  fur  Theoretische  Physik,  TU  Dresden,  D-01062  Dresden,  Germany. 

2  Departement  of  Chemistry,  University  of  Exeter,  Exeter,  EX4  4QD,  UK. 

A  DFT-LDA  Tight  Binding  (TB-DFT)  scheme  for  calculations  of  electronic 
and  geometric  structure  and  properties  of  large  molecules,  clusters  and  solids 
is  presented.  Atomic-like  wave  functions  as  linear  combinations  of  Gaussian 
type  orbitals  are  used  as  basis  functions.  The  effective  Kohn-Sham  (KS)  po¬ 
tential  is  approximated  as  a  superposition  of  KS-potentials  of  neutral  atoms. 
Only  two-centre  matrix  elements  of  the  Hamiltonian  are  considered.  Forces  on 
the  atoms  are  calculated  either  in  a  semiempirical  way  under  consideration  of 
the  electronic  states  or  non  empirically  using  Hellmann-Feynman  forces.  The 
small  computational  effort  of  this  treatment  allows  to  investigate  truly  large 
systems  up  to  size  of  2000  atoms  -  including  molecular  dynamics.  An  inter¬ 
face  allows  to  use  the  results  -  e.g.  electronic  density,  molecular  orbitals,  Fock 
and  overlap  matrices  -  for  property  calculations  using  ’deMon  NMR’  and  ’de- 
Mon  properties’. 

This  method  is  applied  for  the  geometry  optimisation  and  the  calculation 
of  chemical  shifts  of  six  C32  isomers.  The  results  are  compared  with  ’deMon 
KS3p4/NMRlpl’  calculations. 
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Perturbation  theory  of  magnetic  properties 
and  relativistic  corrections  based  on  the 
Levy-Leblond  equation 

Alf  Chr.  Hennum.  Trygve  Helgaker,  and  Wim  Klopper 
Department  of  Chemistry, 

University  of  Oslo, 

P.O.Box  1033  Blindern, 

N-0315  OSLO, 

Norway 

August  10,  1998 


Abstract 

Relativistic  and  magnetic  corrections  are  computed  using  pertur¬ 
bation  theory  based  on  the  Levy-Leblond  equation  —  that  is,  the  non- 
relativistic  limit  of  the  Dirac  equation  —  as  a  zero-order  starting  point. 
The  theory  has  been  applied  through  third-order  to  various  electronic 
states  of  diatomic  one-electron  systems  (Hp  HeH2+,  Hg259+)  and  pre¬ 
liminary  results  for  these  systems  are  presented.  Effects  due  to  the 
finite  size  of  the  nuclei  are  also  studied. 
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Analysis  of  the  Effect  of  Torsion  upon  NMR  Chemical  Shieldings 
and  Spin-Spin  Coupling  Constants  in  Monosaccharide 
Me  P-D-xylopyranoside.  A  Density  Functional  Theory  Approach. 


M.  Hricovmi.  1  O.L.  Malkina, 2  F.  Bizik,^  L.  Turi  Nagy4  and  V.G.  Malkin4 

Institute  of  Chemistry,  1  Computing  Center,2 
Institute  of  Virology f  Institute  of  Inorganic  Chemistry  4 
Slovak  Academy  of  Sciences,  SK-84236  Bratislava,  Slovakia 


Density  functional  theory  based  approaches  were  used  to  calculate  chemical 
shieldings  and  spin-spin  coupling  constants  in  a  monosaccharide  Me  P-D- 
xylopyranoside.  Excellent  agreement  was  found  between  the  computed  and 
experimental  data  for  this  monosaccharide.  Effect  of  torsion  around  the  Cl  - 
01  bond  showed  that  chemical  shifts  of  anomeric  proton,  both  ring  and  01 
oxygens  as  well  as  Cl,  C2  and  the  methyl  carbons  strongly  depend  on  the 
dihedral  angle  where  the  changes  up  to  14  ppm  were  observed.  The  variation  of 
chemical  shielding  of  anomeric  proton  was  found  to  be  0.66  ppm  and  the 
analysis  revealed  that  the  influence  of  01  lone  pairs  upon  its  a  values  play  a 
crucial  role.  The  computed  differences  of  individual  components  of  the 
shielding  tensors  (<Tji)  of  carbons  varied  up  to  50  ppm.  Both  one-bond  and 
three-bond  proton-carbon  coupling  constants  among  anomeric  proton,  anomeric 
and  methyl  carbons  [(1Jci-Hl)  an4  (3JC(Me)-Hl)  ]  respectively,  showed 
strong  dependence  upon  the  torsion  angle  O.  Slight  asymmetry,  as  in  case  of 
^Cl-Hl  couplings,  was  observed  also  in  the  dependence  of  3Jc(Me)-Hl 
couplings  upon  torsion  angle  <£.  The  latter  dependence  was  found  comparable 
with  that  found  experimentally  having,  however,  slightly  larger  values  (~1.5 
Hz)  of  coupling  constants  close  to  the  synperiplanar  conformation. 

Reference 

Hricovini  M.,  Malkina  O.L.,  Bizik  F.,  Turi  Nagy  L.  and  Malkin  V.G. 

"Calculation  of  NMR  Chemical  Shifts  and  Spin-Spin  Coupling  Constants  in 
Monosaccharide  Me  P-D-xylopyranoside  Using  a  Density  Functional  Theory 
Approach",  J.  Phys.  Chem.  A  101(1997)9756-9762. 
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DENSITY  FUNCTIONAL  APPROACH  TO  THE 
STRUCTURES  AND  EPR  PARAMETERS  OF  IMI- 
NOXYL  RADICALS 

A.  R.  Jaszewski.  S.  Berski,  J.  Jezierska 

Faculty  of  Chemistry,  University  of  Wroclaw,  14  Joliot-Curie  St., 
50-383  Wroclaw,  Poland 

Iminoxyls  are  G-type  radicals  characterised  by  a  relatively  large  iso¬ 
tropic  hyperfine  coupling  constant  due  to  a  nucleus  (  a^  =  30  G  )  indi¬ 
cating  a  significant  spin  density  on  nitrogen  in  an  orbital  of  s-containing 
type. 

The  structures  and  EPR  parameters  of  a  set  of  pyridyl  iminoxyl  radi¬ 
cals  are  investigated  by  the  density  functional  (DFT)  approach.  All  the  com¬ 
putations  are  performed  using  Becke’s  three-parametric  hybrid  exchange 
functional  combined  with  the  Lee-Yang-Parr  correlation  functional 
(B3LYP)  and  the  standard  6-31++G(d,p)  basis  set  of  Pople. 

The  molecules  under  study  are  (R)(H)C=N-0  iminoxyls,  where 
R  =  2-pyridyl,  3-pyridyl  and  4-pyridyl.  All  of  them  may  exist  either  in  E  and 
Z  geometry.  Model  considerations  suggest  that  each  isomer  of  2-pyridyl  and 
3-pyridyl  iminoxyl  radicals  has  two  rotameric  forms  (syn  and  anti)  differing 
in  the  orientation  towards  the  C-C  bond  linking  the  >C=N-0*  group  with  a 
pyridyl  ring. 

The  results  of  the  computations  are  in  good  agreement  with  experi¬ 
mental  data  and  offer  a  comprehensive  picture  of  isomeric  effect  on  geomet¬ 
ric  and  electronic  characteristics  (particularly  distinct  in  2-pyridyl  iminoxyl 
radical  through  the  interaction  of  the  pyridyl  nitrogen  with  the  C=N-0*  moi¬ 
ety).  The  calculated  isotropic  and  anisotropic  hyperfine  coupling  constants 
(hfc)  for  and  can  be  used  in  order  to  identify  a  geometric  isomer  of 
iminoxyl  radicals  generated  upon  inclusion  of  2-pyridylaldoxime,  3- 
pyridylaldoxime  and  4-pyridylaldoxime  in  the  channels  of  thermally  acti¬ 
vated  pentasil  zeolite  ZSM-5. 

The  feasibility  of  this  kind  of  computations  for  large  systems  indicate 
that  DFT  methods  seems  to  be  a  useful  tool  for  the  study  of  relationship  be¬ 
tween  structure  and  spectroscopy  properties  of  large  free  radicals. 
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RO  VIBRATION  ALLY  AVERAGED  NUCLEAR  SHIELDING 

CONSTANTS  IN  OCS 


Karol  Jackowski1,  Michal  Jaszuriski2.  Wlodzimierz  Makulski1  and  Juha  Vaara3 


The  nuclear  shielding  constants  in  OCS  are  studied  using  ab  initio  theoretical  methods  and 
gas  phase  NMR  measurements.  The  shielding  surfaces  are  calculated  and  the  rovibrational 
effects  and  the  resulting  temperature  dependence  are  analyzed.  The  temperature  dependence 
of  13C  shielding  in  the  gas  phase  is  determined  experimentally  in  the  range  278 — 373  K.  13C 
is  the  single  nucleus  for  which  the  experimental  data  for  the  temperature  dependence  can 
be  converted  to  a  reference-independent  scale,  and  good  agreement  of  the  measured  and 
calculated  ab  initio  results  is  observed.  For  33S,  we  discuss  a  new,  more  accurate  absolute 
shielding  scale. 
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One-bond  13C-13C  couplings  in  diacetylenes; 
experimental  and  theoretical  studies 

Krystyna  Kamienska-Trela1 ,  Lidia  Kania1,  Wojciech  Schilf , 

and  Irina  Balova2 

I)  Institute  of  Organic  Chemistry,  Polish  Academy  of  Sciences 
Kasprzaka  44,  01-224  Warsaw,  Poland 

2)  Sankt-  Petersburg  University,  Sankt-  Petersburg,  Russia 


Carbon-carbon  spin-spin  coupling  constants  across  one  bond  have  been 
measured  in  several  mono  and  dialkyl  substituted  diacetylenes,  in  dialkyl 
substituted  tetraacetylene  and  in  trimethylsilyl  substituted  diacetylene.  This 
provides  new  examples  of  the  coupling  constants  across  Csp-Csp  single  bonds, 
which,  so  far,  have  been  determined  for  a  very  few  compounds  only.  The 
semi-empirical  INDO  FPT  calculations  have  been  performed  for  mono 
substituted  diacetylenes  with  substituents  varying  along  the  first  row  of  the 
Periodical  Table.  The  calculated  results  revealed  that  one-bond  13C-13C 
coupling  constants  across  C-C  triple  bonds  with  the  substituent  directly  bound 
to  the  bond  concerned  as  well  as  those  across  the  central  single  C-C  bond  are 
significantly  affected  by  substituent  effects. 


P13 


„Normal  vs.  Inverse  Halogen  Dependence44  of  NMR  Chemical  Shifts 
throughout  the  Periodic  Table.  Spin-Orbit  Effects  and  Other  Pecularities 


Martin  Kaupp 

Max-Planck-Institut  fiir  Festkorperforschung,  Heisenbergstr.  1,  D -7 05 69  Stuttgart 
e-mail:  kaupp@vsibml.mpi-stuttgart.mpg.de 

We  have  studied  the  origins  of  „normal  halogen  dependence11  (NHD)  and  „inverse  halogen 
dependence11  (IHD)  of  main  group  and  transition  metal  NMR  chemicals  shifts,  using  DFT-IGLO 
and  DFT-GIAO  calculations  combined  with  a  recently  developed  DFT-based  perturbation 
method  [1]  to  compute  spin-orbit  (SO)  corrections.  The  results  are  discussed  in  the  light  of  a 
recently  delineated  simple  yet  general  qualitative  concept  explaining  such  „SO  chemical 
shifts11. [2]  The  model  is  based  on  the  observation  that  the  SO-induced  spin  polarization  interacts 
with  the  nuclear  magnetic  moments  mainly  via  a  Fermi-contact  (FC)  term,  and  thus  in  part  works 
analogously  to  the  FC  mechanism  of  J-coupling.[2]  Examples  for  this  analogy,  and  for  its  limits, 
are  demonstrated  by  explicit  DFT  calculations  of  the  two  quantities.  NHD  for  main  group  nuclei 
in  the  highest  oxidation  states  is  generally  due  to  SO  effects,  transmitted  via  an  effective  FC 
mechanism,  due  to  the  significant  s-character  of  the  bonding  orbitals.  [2-4]  In  contrast,  no 
efficient  FC  mechanism  is  available  for  p-block  main  group  elements  in  lower  oxidation  states,  as 
their  bonding  is  dominated  by  p-orbitals.  As  a  consequence,  small  SO  effects  and  a  tendency 
towards  IHD  are  found.  Similarly,  small  SO  shifts  and  IHD  are  computed  for  early  transition 
metal  compounds,  due  to  the  predominant  d-character  of  bonding.  NHD  predominates  again  for 
late  transition  metals.  We  find  this  to  be  due  only  in  part  to  SO  effects.  The  main  reason  for  the 
switch  from  IHD  to  NHD  along  the  transition  metal  rows  is  the  different  nature  of  the 
predominant  paramagnetic  shielding  contributions  (involving  ligand-to-metal  charge  transfer 
couplings  for  early,  metal-to-ligand  charge  transfer  couplings  for  late  transition  metals). 

[1]  V.  G.  Malkin,  O.  L.  Malkina,  D.  R.  Salahub  Chem.  Phys.  Lett.  1996,  261,  335. 

[2]  M.  Kaupp,  O.  L.  Malkina,  V.  G.  Malkin,  P.  Pyykko,  Chem.  Eur.  J.  1998,  4,  118. 

[3]  See,  e.g.:  M.  Kaupp,  O.  L.  Malkina,  V.  G.  Malkin  Chem.  Phys.  Lett  1997,  265,  55. 

[4]  See,  e.g.:  M.  Biihl,  M.  Kaupp,  V.  G.  Malkin,  O.  L.  Malkina  J.  Comp.  Chem.,  in  press. 

[5]  M.  Kaupp,  Ch.  Aubauer,  G.  Engelhardt,  T.  M.  Klapotke,  O.  L.  Malkina,  submitted. 
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COMPARISON  OF  AB  INITIO  AND  EXPERI¬ 
MENTAL  ISOTROPIC  HYPERFINE  COUPLING 
CONSTANTS  OF  VARIOUS  SEMIQUINONE 
RADICAL  CATIONS 

A.-M.  Kelterer1,  M.  Sacher1,  D.  Niethammer2 3 4,  and  G.Grampp1 
1  Institut  fur  Physikalische  und  Theoretische  Chemie,  TU  Graz, 
Austria;  2  Institut  fur  Organische  Chemie,  FU  Berlin,  Germany 

A  number  of  semiquinone  radical  cations  of  different  substituted 
N,N,N’,N’-tetraalkyl-p-phenylenediamines  (Fb=Me  [1, 2],  Et  [2, 3, 4],  and 
n-propyl  [3])  were  investigated  experimentally  by  EPR  spectroscopy  and 
theoretically  by  ab  initio  computations. 

Ab  initio  calculations  were  done  by  UHF/3-21G*  and  the  AMI  method 
to  locate  the  low  energy  minima  in  the  potential  energy  surface. 
B3LYP/6-31G*  density  functional  theory  was  used  to  verify  these 
minima  and  to  calculate  the  hyperfine  coupling  constants. 

Comparison  of  experimental  values  with  the  B3LYP/6-3lG*//UHF/ 
3-21G*  calculations  show  a  very  good  correlation  for  the  TMPPD 
radical  cation,  hence  this  level  of  theory  was  used  for  the  other  radical 
cations  too.  Because  of  the  large  number  of  local  minima  with  the 
more  flexible  substituents,  the  ab  initio  data  are  shown  only  for  a  small 
number  of  local  minima. 

[1]  J.  A.  Weil,  J.  R.  Bolton,  J.  E.  Wertz,  Electron  Paramagnetic 
Resonance,  Wiley  &  Sons,  New  York,  1994,  p.  512. 

[2]  G.  Grampp  and  G.  Stiegler,  Z.  Phys.  Chem,  141  (1984)  185. 

[3]  M.  Sacher,  D.  Niethammer,  A.-M.  Kelterer,  and  G.  Grampp,  to  be 
published. 

[4]  J.  P.  Telo,  private  communication. 
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DENSITY  FUNCTIONAL  CALCULATIONS  OF 
HYPERFINE  COUPLING  CONSTANTS  OF  STABLE 
ALANINE  RADICALS 

P.  Lahorte*,  F.  De  Proft**,  B.  Masschaele*,  P.  Cauwels*, 

W.  Langenaeker**,  P.  Geerlings**,  W.  Mondelaers*, 

*Department  of  Subatomic  and  Radiation  Physics, 

University  of  Gent,  Proeftuinstraat  86,  B-9000  Gent,  BELGIUM 
**Eenheid  Algemene  chemie  (ALGC) , 

Free  University  of  Brussels,  Pleinlaan  2,  B-1050  Brussels,  BELGIUM 


The  amino  acid  alanine  is  widely  accepted  as  a  reference  dosimeter  on  the 
basis  of  the  concentration  of  a  stable  radical  which  is  formed  upon 
irradiation  in  the  solid  state  with  high  energy  electrons  or  X-rays.  Recent 
EPR  and  ENDOR  experiments  however,  suggest  the  existence  of  several 
stable  radicals  at  room  temperature. 

In  this  work  Density  Functional  Theory  (DFT)  calculations  (B3LYP 
functional,  variety  of  Pople  basissets)  of  isotropic  hyperfine  coupling 
constants  (HFCC’s)  are  performed  of  both  the  well-known  stable  radical 
and  various  ‘candidate’  radicals  potentially  composing  the  EPR  spectrum  of 
L-alanine.  Calculational  issues  (level  of  geometry  optimisation,  basisset 
quality)  affecting  the  accuracy  of  the  resulting  HFCC’s  are  addressed. 
Although  conclusive  answers  on  the  identity  of  stable  radicals  in  irradiated 
L-alanine  will  require  further  experimental  investigation,  it  is  clear  that  the 
quality  of  the  calculated  coupling  constants  allows  for  differentiation 
between  various  candidate  radicals  of  different  structure  and  geometry. 

The  results  of  this  study  suggest  the  usefulness  of  high  level  standard 
quantum  chemistry  calculations  (DFT)  in  the  determination  of  hyperfine 
coupling  constants  in  molecules  of  real  scientific  and  practical  interest. 
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MCSCF  LINEAR  RESPONSE  CALCULATIONS  OF  CC, 
CF,  AND  FF  SPIN-SPIN  COUPLING  TENSORS 

Perttu  Lantto.  Juha  Vaara,  Jaakko  Kaski,  and  Jukka  Jokisaari 

NMR  Research  Group,  Department  of  Physical  Sciences,  University  of 
Oulu,  P.O.Box  333,  FIN-90571  Oulu,  Finland 


In  NMR  experiments  performed  in  anisotropic  phases,  the  anisotropic  part 
of  the  indirect  KL  spin-spin  coupling  tensor  Ja x  contributes  to  the 
experimental  direct  dipolar  coupling  as  DKlxp  =  DKl  +  ‘/2  Jkl™0-  The 
magnitude  of  Jkl  ms°  should  be  known  if  the  Dklx?  are  used  in  accurate 
determinations  of  molecular  orientation  or  geometry. 

J kl  poses  considerable  requirements  for  basis  sets  and  electron  correlation 
treatment  in  ab  initio  calculations.  We  report  MCSCF  linear  response  [1] 
calculations  for  CC,  CF,  and  FF  couplings  with  all  the  different  physical 
contributions  taken  into  account.  Mostly  large  restricted  active  space  wave 
functions  and  extensive  basis  sets  have  been  used.  The  "JCc  (n  =  1 ,  2,  3) 
couplings  have  been  determined  in  benzene  (C6H6)  [2]  and  in  simple 
hydrocarbons  C2Hm  (m  =  2,  4,  6)  [3],  illustrating  the  dependence  of  the 
coupling  on  the  hybridisation  of  the  carbon  atoms.  The  CF  and  FF  coupling 
tensors  have  been  calculated  in  para- difluorobenzene  (C6H4F2)  [4]  and  in 
the  different  fluoromethanes  CH*F4.*  (k  =  1,  2,  3)  [5]. 

Unlike  in  the  CC  couplings  where  Jr™'™  is  negligible,  it  appears  that  the 
indirect  contribution  to  the  CF  and  FF  couplings  has  to  be  taken  into 
account  in  the  experimental  work. 


[1]  O.  Vahtras,  H.  Agren,  P.  Jprgensen,  H.  J.  Aa.  Jensen,  S.  B.  Padkjajr,  and 
T.  Helgaker,  J.  Chem.  Phys.  96,  6120  (1992). 

[2]  J.  Kaski,  J.  Vaara,  and  J.  Jokisaari,  J.  Am.  Chem.  Soc.  118,  8879  (1996). 

[3]  J.  Kaski,  P.  Lantto,  J.  Vaara,  and  J.  Jokisaari,  J.  Am.  Chem.  Soc.  120, 
3993  (1998). 

[4]  J.  Vaara,  J.  Kaski,  and  J.  Jokisaari,  in  preparation. 

[5]  P.  Lantto,  J.  Kaski,  J.  Vaara,  and  J.  Jokisaari,  in  preparation. 
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Isotope  Effects  on  Structure  and  Dynamics  of  Some 
Radical  Cations  Studied  by  Ab  Initio  and  DFT 
Calculations  and  ESR  Spectroscopy 

S.  Lunell1.  T.  Fangstrom1,  P.  Wang2  and  M.  Shiotani2 
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2 Department  of  Applied  Chemistry,  Faculty  of  Engineering,  Hiroshima 
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Deuterium  isotope  effects  on  conformation  and  ESR  spectra  of  some 
selectively  deuterated  radical  cations  are  studied  using  accurate  ah  initio  and 
density  functional  theory  (DFT)  quantum  chemical  methods,  and  compared 
with  experimental  data  obtained  with  matrix-isolation  ESR  techniques.  The 
temperature  dependence  of  the  ESR  spectra  is  studied  in  the  temperature  range 
4  K  -  100  K.  The  applicability  and  limitations  of  different  computational 
methods  are  illustrated  by  results  on  partially  deuterated  dimethyl  ether  and 
cyclohexane  radical  cations. 


P18 


Dynamical  corrections  to  proton  and  muon  hyperfine  coupling  constants: 
a  comparison  method  for  open-shell  ab  initio  techniques. 

R.  M.  Macrae1,  T.  M.  Briere2 

lMuon  Science  Laboratory,  The  Institute  of  Physical  and  Chemical  Research  (RIKEN), 
2-1  Hirosawa,  Wako-shi,  Saitama  351-01,  Japan 

2  Institute  of  Solid  State  Physics,  University  of  Tokyo, 
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Typical  benchmark  parameters  for  ab  initio  theory  used  in  the  optimisation  of  techniques 
and  basis  sets  are  geometry  (bond  lengths,  bond  angles)  and  harmonic  vibration  fre¬ 
quencies,  and  other  energy  parameters  such  as  heat  of  formation.  For  open-shell  systems, 
another  useful  parameter  is  the  “Fermi  contact”  contribution  to  the  hyperfine  interaction, 
Ax,  proportional  to  the  electronic  spin  density  at  the  nucleus  X  and  normally  computed 
at  the  radical’s  calculated  equilibrium  geometry [?,  ?].  Recent  studies  suggest  that  the 
family  of  density  functional  methods  show  great  promise  for  economical  evaluation  of  Ax 
with  a  precision  comparable  to  that  of  CISD  or  CCSD[?].  However,  dynamical  corrections, 
in  which  the  hyperfine  parameter  is  averaged  over  some  thermally-accessible  segment  of 
the  potential  energy  surface,  greatly  alter  this  picture;  in  certain  cases  such  as  internal  ro¬ 
tation  around  a  partial  double  bond,  density  functional  techniques  show  themselves  to  be 
problematic  [?].  The  present  work  presents  studies  of  torsion  around  a  CO  partial  double 
bond  in  a,/3-unsaturated  ct-hydroxyalkyl  radicals  and  their  isotopomeric  counterparts  in 
which  H  has  been  replaced  by  its  light  isotope,  muonium.  In  general,  the  best  results  (in 
terms  of  comparison  with  experiment)  are  given  by  the  most  computationally-expensive 
techniques,  CISD,  CCSD,  and  MP2.  Additionally,  isotopomeric  studies  of  this  kind  are 
demonstrated  to  be  a  powerful  tool  for  evaluation  of  the  accuracy  of  potential  energy 
surfaces  and  properties  derived  from  the  electronic  wavefunction. 
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The  isotropic  hyperfine  coupling  constants  (hfcc)  for  a  range  of  first-row  transition  metal 
compounds  have  been  investigated  by  means  of  density  functional  theory.  Several  gradient- 
corrected  and  hybrid  functionals  have  been  applied  in  order  to  examine  their  performance  for 
different  systems.  Basis  set  effects  have  also  been  investigated.  A  variety  of  molecules  with 
different  sizes  and  bonding  characters  were  treated,  ranging  from  diatomic  oxides,  nitrides  and 
halides  up  to  penta-  and  hexacoordinated  carbonyl-,  cyano-,  and  nitrosyl-complexes. 

Direct  and  spin-polarization  contributions  to  the  hfcc’s  have  been  analyzed  for  all  complexes. 
Contributions  arising  from  the  spin  polarization  of  2s  and  3s  orbitals  have  been  found  to  increase 
with  the  number  of  unpaired  electrons  in  metal  d-orbitals.  Valence  orbital  contributions  are  not  as 
systematic  and  clearly  depend  on  the  s-character  of  the  SOMO(s),  they  are  positive  for  molecules 
with  no  or  very  weak  s-character  and  negative  otherwise. 

The  effect  of  different  exchange  and  correlation  treatments  on  the  direct  and  indirect 
contributions  to  metal  hfcc’s  has  been  explored.  A  proper  exchange  treatment  has  been  found  to 
be  particulary  important.  Inclusion  of  HF  exchange  in  hybrid  functionals  results  in  more 
pronounced  spin  polarizations  and,  with  a  few  exceptions,  in  weaker  direct  contributions.  As  a 
consequence,  spin  densities  at  the  metal  become  less  positive.  The  general  behaviour  of  different 
functionals  is  as  follows:  B3PW91  gives  metal  hfcc’s  very  close  to  experiment  for  most  of  the 
systems  treated.  Some  high-spin  molecules  with  relatively  low  or  zero  direct  contributions  to  the 
hfcc  are  described  substantially  better  when  50%  HF  exchange  is  included  in  the  density 
functional. 
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NMR  chemical  shifts  in  MNDO  approximation: 
Parameters  and  results  for  H,  C,  N,  and  O. 
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Evaluation  of  the  NMR  chemical  shift  tensor  was  implemented  at  the 
semiempirical  MNDO/d  level  using  gauge-including  atomic  orbitals  and 
analytical  derivative  theory.  Calculations  with  standard  MNDO  parame¬ 
ters  overestimate  the  variation  of  the  paramagnetic  contributions.  Agree¬ 
ment  with  experiment  is  improved  by  reoptimizing  electronic  structure 
parameters  to  reproduce  290  experimental  chemical  shift  values  for  97 
small  molecules  and  ions.  The  reference  set  covers  chemical  shift  ranges 
of  31,  346,  933,  and  1650  ppm  for  H,  C,  N,  and  O,  respectively.  Final 
RMS  deviations  from  experiment  are  0.638,  13.6,  39.6,  and  80.6  ppm 
respectively  for  H,  C,  N,  and  0,  down  from  1.93,  36.8,  120.0,  and  119.0 
ppm  for  original  MNDO.  A  significant  fraction  of  the  total  error  is  due  to 
small  molecules  with  unusual  bonding.  On  a  large  set  of  solution  NMR 
data  for  384  common  organic  molecules,  the  RMS  errors  are  reduced  to 
11.9,  38.8,  and  61.7  ppm  for  C,N,  and  0,  respectively.  Computed  chemi¬ 
cal  shifts  are  sensitive  to  the  molecular  geometry.  Bond  lengths  accurate 
to  at  least  0.05 A  are  normally  required  for  reasonable  agreement  with 
experiment.  Three-center  terms  typically  contribute  a  few  ppm  to  the 
total  chemical  shift  for  all  four  elements  studied  and  are  essential  for  a 
qualitatively  correct  description  of  XH  chemical  shifts  and  NICS.  These 
terms  are  less  important  for  C,  N,  and  0  and  can  usually  be  omitted  with¬ 
out  significantly  degrading  the  results.  Notable  exceptions  are  molecules 
with  cumulated  double  bonds  and  systems  with  long-range  effects  due 
to  aromatic  rings.  The  overall  quality  of  the  results  obtained  with  the 
NMR-specific  parametrization  is  competitive  with  non-correlated  ah  ini¬ 
tio  calculations. 


P21 


COMPARATIVE  STUDIES  OF  THE  EXPERIMENTAL  AND  CALCULATED 
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The  subject  of  this  paper  was  analysis  of  the  electron  density  and  populations  of  the  s  and  p 
orbitals  in  inorganic  compounds  as  well  as  complexes  containing  such  atoms  as  Sn,  Sb,  Al,  I 
and  Cl  with  different  organic  ligands.  The  use  of  NQR  and  Mossbauer  spectroscopies  enabled 
simultaneous  analysis  of  changes  in  the  s-and  p-orbital  electron  density  occurring  upon  complex 
formation.  Results  were  referred  to  the  distributions  of  electron  density  calculated  by  ab  initio 
methods.  The  calculations  were  carried  out  with  the  program  packages  GAMES S,  MOLCAS, 
HONDO  and  GAUSSIAN  94,  run  on  an  IBM  RS/6000  work  station  in  the  double  polarised 
TZ2P  basis  set.  Results  were  obtained  on  the  HF  4-3 1G,  6-3 11G**  and  ccpVZ  levels  of  the 
theory.  The  quadrupole  coupling  constants  for  Sn,  Sb,  Al,  I,  Cl  and  Br  were  calculated  on  the 
basis  of  the  eigenvalues  of  the  electric  field  gradient  tensor  using  Townes-Dailey 
approximation. 

The  NQCC  constants  for  the  central  atoms  (Sn,  Sb  and  I)  calculated  by  the  program  packages 
HONDO  and  GAUSSIAN  using  the  extended  basis  6-311  G**  with  pseudo-potential  were 
much  lower  than  the  experimental  data.  Calculations  performed  by  MOLCAS  and  GAMESS 
packages  using  the  4-3 1G  basis  set  led  to  NQCC  much  closer  to  the  experimental  ones.  The 
frequencies  calculated  for  the  terminal  halogen  atoms  are  well  correlated  with  the  experimental 
results  in  6-3 1 1G**  basis.  Also,  a  good  agreement  with  the  experiment  was  obtained  for  NQCC 
values  calculated  for  the  central  atoms  of  some  compounds  by  the  PM3  and  AMI  methods 
using  Townes-Dailey  approximation.  We  concluded  that  application  of  the  pseudo-potential  by 
HONDO  and  GAUSSIAN  packages  did  not  lead  to  appropriate  NQCC  values  for  the  non¬ 
transition  atoms  containing  many  electrons  such  as  Sn,  Sb  and  I.  Analysis  of  the  chemical  shifts 
in  the  Mossbauer  spectra  of  the  model  complexes  (in  which  the  bond  lengths  between  central 
atom  and  organic  ligand  were  changed)  containing  either  antimony  or  tin  led  to  the  conclusion 
that  in  the  case  of  Sb  atoms  the  shifts  were  caused  by  a  change  in  the  5  s  -orbital  population 
while  the  shifts  attributed  to  cis-SnClff^  complexes  could  be  related  to  the  shielding  of  5p- 
orbitals.  The  different  behaviour  of  Sn  and  Sb  is  a  consequence  of  a  different  hypervalency  of 
bonds  determined  by  a  greater  contribution  of  electron  density  coming  from  the  s-orbitals  of  the 
central  atom  in  the  chemical  bond  between  the  ligand  and  Sn  than  between  the  ligand  and  Sb. 


The  Electric  Field  Dependence  Of  Nuclear  Spin-Spin  Coupling 


M.  Grayson  and  W.T.  Raynes 


Dept,  of  Chemistry,  The  University  of  Sheffield,  Sheffield,  S3  7HF,  UK 


Ab  initio  calculations  of  the  electric  field  dependences  of  nuclear  spin-spin 
coupling  have  been  carried  out  for  the  first  time.  The  systems  under  study 
include  each  of  several  small  molecules  placed  in  a  uniform  electric  field.  The 
derivatives  of  the  coupling  up  to  second  order  in  the  electric  field  have  been 
determined.  The  number  of  independent  coefficients  depends  on  the  "site 
symmetry"  of  the  coupled  nuclei  in  relation  to  the  molecule  as  a  whole.  For  the 
case  of  J(H,H)  in  molecules  such  as  water  and  methanal  allowance  must  be  made 
for  the  exchange  symmetry  of  the  proton  pair. 

Results  will  be  presented  for  the  water  molecule  and,  it  is  hoped,  a  number  of 
other  small  molecules.  For  water  calculations  at  the  SCF  and  MCSCF  levels 
show,  as  expected,  that  electron  correlation  effects  are  very  important.  For 
example,  for  the  proton-proton  coupling  dJ(iso)/dFz  is  -5.45  Hz  (a.u.)'1  at  the  SCF 
level  and  +  34.15  Hz  (a.u.)'1  at  the  MCSCF  level  where  the  z-axis  is  the  C2  axis. 
Almost  all  of  this  change  is  due  to  the  effect  of  electron  correlation  on  the 
derivative  of  the  Fermi  contact  term. 
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Magnetic  Field  Density  Functional  Theory 

Freddie  R.  Salsburv.  Jr.tt.  and  Robert  A.  Harris^ 
tMaterials  Sciences  Division,  Lawrewnce  Berkeley  National  Laboratory 
^Department  of  Chemistry,  University  of  California,  Berkeley 

Over  the  past  2  years  an  alternative  form  of  density  functional  theory,  Magnetic 
Field  Density  Functional  Theory  (BDFT),  has  been  formulated  (1).  This  theory,  in 
principle,  can  allow  for  simpler  computation  and  interpretation  of  NMR  parameters,  though 
much  computation  work  remains  to  be  performed.  The  basic  elements  of  this  theory  will 
be  presented. 

(1)  F.  R.  Salsbury,  and  R.  A.  Harris,  J.  Chem.  Phys.,  108  (15)  6102  (1998)  and 
references  therein. 
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Determination  of  Proton  Hyperfine  Tensors  in  Single  Crystals  of 
D.  vulgaris  Miyazaki  F  [NiFe]  Hydrogenase.  Pulsed-ENDOR 
Spectroscopy  and  DFT  Calculations. 

M.  Stein,  0.  Trofanchuk,  M.  Brecht,  F.  Lendzian,  Y.  Higuchi§ , 

W.  Lubitz 

Max-Volmer-Institut  fiir  Biophysikalische  Chemie  und  Biochemie,  Technische  Universitat 
Berlin,  Strasse  des  17.  Juni  135,  D-10623  Berlin,  Germany 
§  Division  of  Chemistry,  Graduate  School  of  Science,  Kyoto  University,  Sakyo,  Kyoto 

606-01,  Japan 

By  Electron  Nuclear  Double  Resonance  (ENDOR)  spectroscopy  the  interac¬ 
tions  between  electron  spin  and  (ligand)  nuclear  spins  in  paramagnetic  metal 
complexes  can  be  determined.  The  investigation  of  protein  single  crystals 
allows  the  measurement  of  complete  hyperfine  tensors,  i.e.  isotropic  and 
anisotropic  contribution,  of  the  active  center  of  paramagnetic  biomolecules. 
The  orientation  of  the  g-tensor  in  the  NiFe  cluster  of  Miyazaki  F  has  already 
been  determined  [1].  The  ’’ready”  Ni-B  state  of  the  [NiFe]  hydrogenase  of 
D.  vulgaris  Miyazaki  F  [2]  was  investigated  in  single  crystals  by  X-Band 
(9.5  GHz)  1H  pulsed-ENDOR  spectroscopy.  We  were  able  to  determine  the 
complete  hyperfine  tensors  of  three  protons  in  the  vicinity  of  the  active  cen¬ 
ter.  The  analysis  of  the  experimental  data  gave  two  proton  hyperfine  tensors 
with  large  isotropic  contributions  (11  and  10  MHz)  and  one  smaller  almost 
purely  dipolar  hyperfine  interaction  ( a;so  =  5  MHz).  Based  on  comparison 
of  the  anisotropic  dipolar  contributions  of  the  respective  hyperfine  tensors 
with  the  X-ray  structure  the  two  /T-cysteinyl  protons  of  the  bridging  cysteine 
Cys549  are  assigned  to  the  larger  two  hyperfine  interactions  while  the  third 
mainly  dipolar  interaction  is  assigned  to  a  protonated  bridging  ligand  X,  i.e. 
X=OId-  or  SH“.  DFT  calculations  on  a  large  cluster  model  of  the  active 
center  support  this  assignment.  The  spin  density  is  in  agreement  with  ex¬ 
perimental  data  when  the  bridging  ligand  is  assumed  to  be  OH-. 
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Exchange  broadening  dependence  of  EPR  spectra  on 

NITROXIDES’  ELECTRONIC  STRUCTURE  AND  MEDIA  PROPERTIES 


Janez  Strancar,  Marjana  Nemec,  Milan  Schara 
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Nitroxide  radicals  are  widely  used  spin  probes  exploring  the  complexity  of 
biological  systems,  since  their  electronic  structure  and  therefore  EPR  spectra 
are  quite  sensitive  to  the  local  environmental  properties. 

In  this  context  we  studied  the  EPR  spectra  lineshapes  influenced  by  the  spin 
exchange.  This  is  one  of  the  effect  arising  from  electronic  structure 
variation.  We  varied  the  type  of  residue  attached  on  oxyl-piperidin  ring, 
their  concentration  and  the  solution  media.  Three  nitroxide  labels  were 
inspected  -  tempo,  tempone  and  glusl  in  different  solutions  -  water,  ethanol 
and  acetone,  chosen  to  cover  the  whole  interval  of  hyperfme  tensor  changes. 
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With  concentration  dependence  we  measured  the  exchange  frequency  which 
is  linearly  dependent  on  the  spatial  integral  of  overlapping  orbitals.  This 
overlap  integral  differs  by  radical  type  and  media  polarity  due  to 
delocalization  of  unpaired  n  electron  over  piperidine  ring. 

Temperature  dependence  was  measured  to  check  the  effect  of  dipolar 
broadening  which  in  case  of  small  nitroxide  molecules  and  therefore  fast 
rotational  motions  usually  vanishes. 

We  also  measured  the  spectra  in  the  presence  of  oxygen  and  without  it  to 
exclude  the  oxygen  broadening  effect,  what  is  especially  important  in 
nonpolar  media  due  to  higher  oxygen  solubility. 
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The  Aromatic  Pathway  of  Chlorins  and 
Porphyrins:  A  Computational  Study 

Jonas  Juselius  and  Page  Sundholm 

Department  of  Chemistry,  P.O.Box  55,  FIN-00014  Helsinki,  Finland 

Abstract 

The  aromatic  pathway  in  chlorins  and  porphyrins  has  been  studied  by 
performing  Hartree-Fock  (SCF)  calculations  of  the  nuclear  magnetic  shiel¬ 
ding  constants  using  TURBOMOLE  [1-2].  The  accuracy  of  the  obtained 
SCF  shieldings  was  checked  at  second-order  Mpller-Plesset  level.  Our 
calculations  show  that  the  preferred  aromatic  pathway  does  not  include 
the  NH  bridge  as  recently  suggested  by  Cyrahski  et  al.  [4];  instead  it 
forms  a  classical  [18]annulene  structure  with  the  NH  groups  as  relatively 
inert  bridging  units. 
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SPIN-ORBIT  EFFECTS  ON  NMR  SHIELDING  AND  SPIN- 
SPIN  COUPLING  TENSORS  BY  QUADRATIC 
RESPONSE  THEORY 

J.  Vaara,1 2 3 4 5 6  K.  Ruud,b  O.  Vahtras,c  B.  Minaev, c  H.  Agren,c  and  J.  Jokisaari3 

aNMR  Research  Group,  Department  of  Physical  Sciences,  University  of 
Oulu,  P.O.  Box  333,  FIN-9057 1  Oulu,  Finland 
bDepartment  of  Chemistry,  University  of  Oslo,  Blindern,  P.O.  Box  1033, 

N-0315  Oslo,  Norway 

department  of  Physics  and  Measurement  Technology,  Linkoping 
University,  S-58  183  Linkoping,  Sweden 

Due  to  their  importance  for  the  NMR  of  systems  containing  heavy  elements, 
the  calculation  of  relativistic  effects  -  the  spin-orbit  (SO)  effect  in  particular 
-  on  NMR  parameters  has  gained  significant  interest  lately.  The  SO-effect 
causes  the  heavy-atom  shift  of  the  NMR  signals  of  light  nuclei  located  in 
the  neighbourhood  of  heavy  atoms  [1-3].  Ab  initio  methods  based  on  SCF 
and  MCSCF  reference  states  and  the  quadratic  response  theory  [4]  for 
calculating  the  S O-corrections  to  shielding  and  coupling  tensors  have 
recently  been  developed  [5,6].  We  discuss  the  theory  in  detail  and  give 
explicit  expressions  for  the  tensor  elements.  Numerical  examples  will  be 
provided  for  the  HX  and  CH3X  (X  =  F,  Cl,  Br,  I),  H2X  (X  =  O,  S,  Se,  Te), 
and  XeF2  systems,  as  well  as  the  bond  length  dependence  of  the  SO- 
contribution  to  *H  shielding  in  HI. 
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DICHROISM  IN  THE  ELECTROMAGNETIC 
RESPONSE  OF  SUPERCONDUCTORS 
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Institut  fiir  Theoretische  Physik,  Universitat  Wurzburg, 
Am  Hubland,  D-97074  Wurzburg,  Germany 


A  relativistic  generalization  of  the  Bogolubov-deGennes  equations  will  be 
presented  where  the  particle  and  hole  amplitudes  are  Dirac  spinors.  In  the 
weakly  relativistic  limit  one  obtains,  besides  the  usual  spin-orbit,  Darwin  and 
kinetic  energy  corrections,  an  additional  “spin-orbit”  and  an  additional  “Dar¬ 
win”  term  [1],  The  latter  terms  are  present  in  superconductors  only  and  in¬ 
volve  the  pairing  field  in  place  of  the  electrostatic  potential.  In  view  of  the  fact 
that  many  unconventional  superconductors  have  at  least  some  heavy  elements 
in  the  lattice,  the  significance  of  these  new  terms  for  possible  superconduct¬ 
ing  mechanisms  in  the  heavy-fermion  and  the  high-Tc  compounds  will  be 
discussed. 

As  a  first  application,  a  microsopic  theory  of  dichroism  in  the  magneto-optical 
properties  of  superconductors  will  be  presented  [2,3].  Four  distinct  mech¬ 
anisms  for  dichroism  in  superconductors  are  identified.  Two  are  modifica¬ 
tions  of  mechanisms  known  from  the  normal  state,  and  two  are  novel  effects 
found  in  superconductors  only.  The  theory  is  illustrated  by  numerical  calcula¬ 
tions  for  a  simple  model  system.  The  interplay  between  relativistic  symmetry 
breaking  and  superconducting  coherence  is  found  to  give  rise  to  a  variety  of 
new  effects,  not  known  from  dichroism  in  the  normal  state. 

[1]  K.  Capelle,  E.K.U.  Gross,  Phys.  Lett.  A  198  (1995),  261. 

[2]  K.  Capelle,  E.K.U.  Gross,  B.L.  Gyorffy,  Phys.  Rev.  Lett.  78  (1997),  3753. 

[3]  K.  Capelle,  E.K.U.  Gross,  B.L.  Gyorffy,  Phys.  Rev.  B  58  (1998),  473. 
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Density  Functional  Calculations  on  Sodium  Doped 

Sodium  Sodalite 

Rene  Windiks.  Joachim  Sauer 

Humboldt-Universitat  zu  Berlin,  Institut  fur  Chemie,  Arbeitsgruppe 

Quantenchemie, 

Jagerstr.  10/11,  D-10117  Berlin,  Germany 

When  dehydrated  sodium  aluminosilicate  sodalite,  (Na+)3(AlSi07)3, 
is  doped  with  sodium  atoms,  paramagnetic  (Na+)4e_=  Na*3  clusters 
are  formed  in  its  cages.  Due  to  the  isotropic  hyperfine  (Fermi  contact) 
interactions  the  23Na,  27 A1  and  29Si  MAS  NMR  spectra  show  several 
paramagnetically  shifted  resonances.  The  corresponding  spin  densities 
at  the  Si  and  A1  nuclei  of  the  framework  and  at  the  Na  nuclei  of  the 
diamagnetic  (Na+)3  clusters  are  in  the  range  of  10-5  and  10-4  atomic 
units,  while  at  the  Na  nuclei  of  the  paramagnetic  Na|3  clusters  the  spin 
densities  are  of  the  order  of  10_1  atomic  units. 

We  perform  density  functional  calculations  to  examine  the  assump¬ 
tion  that  the  observed  multiple  lines  in  the  27A1  and  29Si  NMR  spectra 
are  due  to  the  fact  that  the  spin  density  at  a  specific  T  site  (T  =  Al, 
Si)  depends  on  the  number  of  surrounding  paramagnetic  Na|3  clusters 
[1].  We  also  rationalize  the  presence  of  spin  density  at  the  nuclei  of  the 
diamagnetic  (Na+)3  clusters. 

Calculations  are  performed  on  periodic  structures  as  well  as  on  cluster 
models.  The  periodic  calculations  include  as  many  as  88  atoms  (2x1x1 
cells)  and  employ  the  DSolid  and  WIEN95  codes.  For  the  calculations  on 
cluster  models  it  is  important  that  (i)  the  composition  of  the  models  cut 
out  of  the  periodic  structure  is  similar  to  the  stoichiometric  composition 
of  the  crystal,  (ii)  the  model  is  made  electrically  neutral  by  adding  an 
array  of  point  charges  and  (iii)  the  symmetry  of  the  environment  of  the 
magnetic  nucleus  in  the  model  is  the  same  as  in  the  crystal. 

[1]  Engelhardt,  G.;  Feuerstein,  M.;  Sieger,  P.;  Markgraber,  D.;  Stucky, 

G.  Srdanov,  V.  J.  Chem.  Soc.  Chem.  Commun.  6,  729  (1996). 
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